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^ : ABSTRACT 

We present the first high spatial resolution X-ray study of NGC 2244, the 2 Myr old 
,^ ] stellar cluster immersed in the Rosette Nebula, using the Chandra X-ray Observatory. 

»j^| Over 900 X-ray sources are detected; 77% have optical or FLAMINGOS near-infrared 

P ' (NIR) stellar counterparts and are mostly previously uncatalogued young stellar cluster 

members. All known OB stars with spectral type earlier than Bl are detected and the 
j^ I X-ray selected stellar population is estimated to be nearly complete between 0.5 and 

3 M0. The X-ray luminosity function (XLF) ranges from 29.4 < logL^ < 32.0 ergs s~^ 
ps| . in the hard (2 — 8 keV) band. By comparing the NGC 2244 and Orion Nebula Cluster 

^ I XLFs, we estimate a total population of ~ 2000 stars in NGC 2244. A number of further 

^1 ' results emerge from our analysis: 1) The XLF and the associated X-band luminosity 

f^ . function indicate a normal Salpeter initial mass function (IMF) for NGC 2244. This is 

^N I inconsistent with the top-heavy IMF reported from earlier optical studies that lacked 

a good census of < 4M0 stars. 2) The spatial distribution of X-ray stars is strongly 
concentrated around the central 05 star, HD 46150. The other early O star, HD 46223, 
("^ [ has few companions. The cluster's stellar radial density profile shows two distinctive 

structures: a power-law cusp around HD 46150 that extends to ~ 0.7 pc, surrounded 
by an isothermal sphere extending out to 4 pc with core radius 1.2 pc. This double 
structure, combined with the absence of mass segregation, indicates that this 2 Myr old 
d I cluster is not in dynamical equilibrium. Our results will strongly constrain models of 

the cluster formation process. The spatial distribution of X-ray selected K-excess disk 
stars and embedded stars is asymmetric with an apparent deficit towards the north. 3) 
The fraction of X-ray-selected cluster members with A'-band excesses caused by inner 
protoplanetary disks is 6%, slightly lower than the 10% disk fraction estimated from the 
FLAMINGOS study based on the NIR-selected sample. This is due to the high efficiency 
of X-ray surveys in locating disk- free weak- lined T Tauri stars. 4) X-ray luminosities 
for 24 stars earlier than B4 confirm the long-standing log{Lx/Lhoi) ~ —7 relation. The 
Rosette OB X-ray spectra are soft and consistent with the standard model of small- 
scale shocks in the inner wind of a single massive star. 5) About 50 intermediate-mass 
(2 < M < 8 Mq) cluster members are identified; they exhibit a wide range of X-ray 
luminosities consistent with previously studied samples of Herbig Ae/Be stars. Their 
A'-excess disk fraction is ~10%, indicating that the Herbig Ae/Be phenomenon is rare 
or short-lived. 

Subject headings: Open clusters and associations: individual (NGC 2244) -ISM: individual 
(Rosette Nebula) - stars: formation - stars: mass function - stars: pre-main sequence - X-Rays: 
stars 
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1. Introduction 

During their evolution from Class I pro- 
tostars to zero-age-main-sequence (ZAMS) 
stars, young stellar objects are readily iden- 
tified in X-rays due to their highly ele- 
vated X-ray emission compared to the older 
Galactic stellar population (see reviews by 
Feigelson k. Montmerlelll999l : lFavata fc Micela 



20031 : iFeigelson et al.ll2007l '). High spatial res- 
olution Chandra observations of well-known 
Galactic star forming regions (e.g., the Chan- 
dra Orion Ultradeep Project, hereafter COUP, 
Getman et al. 2005; RCW 38, Wolk et al. 
2006; Cepheus B, Getman et al. 2006), more 
distant molecular cloud and HII region com- 
plexes (e.g., NGC 6334, Ezoe et al. 2006; NGC 
6357, J. Wang et al. 2007; M 16, Linsky et al. 
2007; M 17, Broos et al. 2007), the Galactic 
Cent er (e.g., the Arches and Quintuplet clus - 
ters, lO. Wang et"all l2006l : IMuuo et al.l I2OO6I I , 
and ex tragalactic star formin g regions (e.g., 
30 Dor. lTownsley et al.ll2006al )bl). have greatly 
advanced our knowledge of star formation 
processes in these regions. 

Moreover, these studies demonstrate the 
unique power of studying star formation in 
the X-ray band. Besides the high energy 
phenomena, other new information about the 
young clusters such as population, member- 
ship, and star formation environs can be ob- 
tained when combining X-ray detections with 
knowledge obtained from longer wavelength 
studies. For example, although X-ray-selected 
samples contain a very small fraction of ex- 
tragalactic sources and Galactic field stars, 
X-ray sources that have optical and infrared 
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(IR) counterparts are mostly cluster members. 
This is in contrast to the high percentage of 
non-members in the optical/IR images where, 
except for those sources with massive dusty 
disks, membership for individual stars gener- 
ally has to be ascertained via spectroscopy. 
In initial mass function (IMF) studies of more 
distant, high mass star forming regions, tradi- 
tional optical measurements are significantly 
encumbered by large reddening and member- 
ship confusion. In comparison to the mas- 
sive stars, the lower mass pre-main-sequence 
(PMS) stellar populations are much less ac- 
cessible and thus much harder to evaluate. In 
the past. Ha emission was used in general as a 
youth and membership indicator. However, in 
addition to the observational challenges due to 
prevalent bright H^ nebulosity in HII regions, 
the emission itself requires accretion activ 



ity in prot oplanetary disks ([Muzerolle et al 



19981 . I2OOII ) and thus is susceptible to the disk 



evolutionary stages. 

It has long been recognized that X-ray 
emission circumvents these problems and is 
very effective f or securing PMS membership o f 
young clusters (IFeigelson &: Montmerlelll999l ) . 
Modern X-ray observatories like Chandra and 
XMM-Newton enable identification of hun- 
dreds of individual members from their X-ray 
emission. In this work and subsequent papers, 
we report Chandra studies of a well-known 
star forming complex, concentrating here on 
a new census of the low mass cluster members 
and new knowledge of the IMF arising from 
the X-ray perspective. 

The Rosette star-forming complex, situ- 
ated in a large star formation site in the 
Perseus spiral arm, provides an ideal testbed 
for studying sequential formation of clusters 
due to the favorable orientation of its mor- 
phological components, consisting of an ex- 
panding blister HII region on the edge of 
a giant molecular cloud oriented perpendic- 
ular to the line-of-sight. Thorough review 



of past and present research on this popular 
an d important star forin ation region appear 



m 



Townslev et alJ (120031^ fh ereafter TFM03) 



and lRoman-Zuhiga &: Lada I ( 20071 ) (hereafter 
RL07). Here we highhght what is most rele- 
vant to this work. 

The Rosette Nebula (= Sharpless 275 = W 
16 = NGC 2237-2239, NGC 2244, and NGC 
2246) is a large HII region at the tip of the 
Rosette Molecular Cloud (RMC). In both ra- 
dio and optical images, it shows a prominent 
ring- like morphology, with a cluster of i onizing 
young stars located iii the c entral hole (|Celnik 



19851 : iTownslev et al. 
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J). The large-scale 
^^_ ) and CO emission 

map (|Hever et al...2006 ) clearly show a similar 
annular morphology extending into the molec- 
ular cloud (Figure [1^). 

Extending to the southeast of the Rosette 
Nebula, the RMC is an elongated giant molec- 
ular cloud with ^ lOfM© of gas and dust 
(JBlitz fc ThaddeusI ligsd ). Multiwavelength 
observations in mid-IR and r adio show clumpy 
struc ture in the RMC (jWilliams fc Blitz 
19981 ). Embedded star clusters have been 
revealed in the denses t parts through near- 



19971 : 



IR imaging surveys (iPhelps fc Ladal 
Roman-Zufiiga et al.l l2005l . l2007l ). F igure [Tb 
demon strates the association between lPhelps fc 
(jl997l ) IR clusters and the mol ecular clumps 
where the CO emission peaks. iLada fc Lada 
( 20031 ) have shown that embedded clusters are 
physically associated with the most massive 
and dense cores in molecular clouds, based 
on systema tic and coordinated surveys (e.g., 
Ladalll992l ^. 

In the X-ray band, previous imaging study 
of the Rosette Complex was hampered by 
low spatial resolution. An early Einstein 
Observatory study detected a few individ- 
ual O stars and extended ~ 2 keV X-ray 
emiss io n at the center of the nebula (JLeahv 
igsi). iBerghofer fc ChristianI (|2002l ^ consid- 



emitting low mass PMS stars in the cluster 
and concluded that apparently diffuse emis- 
sion seen by ROSAT coul d be explai n ed by 
unresolved point sources. IChen et al.l ( 20041 ) 
analyzed the same ROSAT data set and at- 
tributed the brightest X-ray sources to mas- 
sive stars, act ive T Tauri stars, and fore - 
ground stars. ICregorio-Hetem et al.l ( 19981 ) . 
studying a ROSAT observation of the RMC, 
reported faint X-ray point sources associ- 
ated with T Tauri stars and Herbig Ae/Be 
(HAeBe) stars and X-ray "hot spots" from 
unresolved embedded low mass star clusters. 
In TFM03, the first high resolution Chandra 
X-ray image mosaic of this high-mass star 
forming region was presented and soft diffuse 
X-ray plasma {kT ~ 0.06 and ~ 0.8 keV) 
with luminosity L^,. — 6 x 10^^ ergs s~^ was 
detected in the HII region. It was attributed 
to a combination of fast O star winds and 
unresolved T Tauri stars. 

The young stars powering the Rosette HII 
region are members of the massive open clus- 
ter NGC 2244. Figure lb shows a Digital Sky 
Survey image (59' x 59') of this region with 
a few bright stars and interesting objects la- 
beled. Despite its apparent low concentra- 
tion of stars in the optical, NGC 2244 con- 
Lada tains ~ 30 early-type stars between 04V and 
B3V (Table 6 in TFM03), whose cluster mem- 
berships are secured from deep photometric 
study along with proper mot i on data and 
spect roscopy (jVerschuerenlll99ll : lPark &: Sung 



2OO2I ). No other known massive young stel- 
lar cluster within 2 kpc, other than RCW 



38 (IWolk et al.ll2006l l and M 17 (Broos et al. 
2007), is comparably rich. The PMS mem- 
bers were largely unknown from optical stud- 
ies; only a handful of Ha emission objects 
have been identified as young PMS mern - 
bers dPark fc Sung|l2002l: IlI &: Rectoil l2004l) . 



ered the integrated contribution from X-ray 



Berghofer fc ChristianI (|2002l ) presented 138 
sources selected from ROSAT VSPC and HRI 
detections. Although incomplete, they re- 



vealed previously unknown PMS stars to K 
spectral types. 

The distance to NGC 2244 has been mea- 
sured in many visual photometry studies 
and ranges between 1 .4 kpc and 1.7 kpc. 
Hensberge et al.l (2000) derive a distance of 
1.39 lb 0.1 kpc to an eclipsing binary mem- 
ber V578 Mon, using a novel Fourier spectral 
disentangling technique. The n iain sequence 
(MS) turn-off age estimated by iPark &: Sung 
(2002) is 1.9 Myr, which is consistent with 
the in ferred age from y578 Mon, 2.3 it 0.2 



Myr ( Hensberge et al.l l200d ). This makes 



NGC 2244 the youngest cl uster within the 



large r Mon 0B2 association ([Hensberge et al 
2000). We adopt a distance of 1.4 kpc 



and a cluster age for NGC 2244 of 2 Myr 
throughout our studies; these are consis- 
tent with TFM03. Note that a larger dis- 
tance value d = l.G kpc has been used re- 
centl y by other researchers (IPark &: Sung 
2OO2I : iRoman-Zuniga et al.l l2007l l. This dis- 
crepancy in distance affects our derived X-ray 
luminosities (logL^;) by only 0.1 dex. At the 
distance of 1.4 kpc, 1' corresponds to 0.4 pc. 

Whereas the TFM03 paper was dedicated 
to study of the diffuse X-ray emission in the 
Rosette HII region, we present here an X- 
ray point source study of the Rosette complex 
based on TFM03 data and a new 75 ks Chan- 
dra observation centered on NGC 2244. The 
field of view of our mosaic Chandra fields are 
outlined by the polygons in Figure [1^. We 
separate our study into a series of four pa- 
pers with different astrophysical emphasis. In 
this work (Paper I), we report Chandra ob- 
servations of the NGC 2244 cluster and the 
Rosette HII region and study the young stel- 
lar population in detail. In an upcoming pa- 
per (Wang et al., in preparation; Paper II), we 
describe Chandra observations of the embed- 
ded clusters in the RMC, aiming to investi- 
gate cluster formation in a sequential manner 
and to test whether molecular clumps pref- 



erentially forming embedded clusters of low- 
mass stars make up the fundamental building 
blocks of star formation in molecular clouds. 
The westernmost Chandra field, designed to 
study triggered star formation and the X-ra; 
detection of a twin cluster to NGC 2244 



20051 : iRoman-Zuhiga et al.ll2007l ). will be pre- 
sented in Paper III (Wang et al., in prepara- 
tion) . A detailed analysis of the diffuse X-ray 
emission in the HII region will appear in Pa- 
per IV (Townsley et al., in preparation). 

New mid-IR observations are also con- 
tributing to our knowledge of the Rosette 
complex. A shallow Spitzer Space Tele- 
scope survey of NG C 2244 has been reported 
( Balog et al.ll2007l ) and a deep Spitzer survey 
for disk emission from Chandra stars in NGC 
2244 is also underway (PI: Bouwman). To- 
gether with the Spitzer MIPS coverage of the 
massive cores in the RMC (PI: Bonnell) and 
a new Spitzer program (PI: Rieke), virtually 
the entire nebula and the molecular cloud will 
be completely mapped. 

This paper is organized as follows. First, 
we describe the Chandra observations and 
data reduction in §2. In §3, we identify 
the X-ray sources with optical and infrared 
counterparts, and evaluate the fraction of 
contaminants through simulations utilizing 
stellar population synthesis model and the 
logA'^-logS' distribution for extragalactic X- 
ray sources. §4 is devoted to global properties 
of the NGC 2244 cluster such as the X-ray 
luminosity function, the initial mass function, 
the i^-band luminosity function, spatial struc- 
tures, mass segregation, and the X-band ex- 
cess disk fraction among the X-ray detected 
stars. We present collective properties of in- 
teresting X-ray sources in §5, ending with a 
summary in §6. 



2. Chandra Observations and Data Re- 
duction 

The Rosette complex was observed with 
the Imaging Array of the Chandra Advanced 
CCD Imaging Spectrometer (ACIS-I). The 
ACIS-I field of view is 17' x 17' in a single 
pointing and a mosaic observation was de- 
signed to best image the ionizing cluster, cap- 
ture the interface between the photoionized 
gas and the cold neutral material, step into 
the dense molecular cloud, and study its re- 
cently reported secondary cluster. As shown 
in Table 1, the entire observation consisted of 
four ~20 ks ACIS-I snapshots in January 2001 
(TFM03, Figure 2), a deep 75 ks ACIS-I im- 
age in January 2004 centered on the 05 star 
HD 46150 in NGC 2244 (Figure [3i), and one 
20 ks ACIS -I pointi ng at the twin cluster to 
NGC 2244 (JLil koO,^ ) in 2007. The image mo- 
saic covers a ~ 1° x 0°.25 field of the Rosette 
Nebula and RMC. All images were taken in 
standard "Timed Event, Faint" mode with 3 
pixel X 3 pixel event islands except ObsID 
3750 and ObsID 8454, which used the "Very 
Faint" mode (5 pixel x 5 pixel event islands). 

We follow the same customized data re- 
duction and source extraction described in 



TFM03. 1 Wang et al.1 (|2007l ^. and iBroos et al. 



( 20071 ). The processing of Level 1 data is pre- 
sented in detail in Appendix B of TFM03; 
the same reduced dataset used in that study 
(Rosette Field 1-4), augmented by the deep 
observation (Rosette Nebula/NGC 2244) re- 
duced following TFM03, were used here for 
further analysis. With slightly different roll 
angles, we reprojected the ObsID 1874 data 
and merged them with the ObsID 3750 field. 
Figure [2^ shows the merged 94 ks ACIS im- 
age of NGC 2244 overlaid with source extrac- 
tion regions (see details below). Many point 
sources are visible; this is further illustrated 
in the smoothed X-ray composite image (Fig- 
ure [2)3) for the merged fields created with the 



CIAO tool csmooth (JEbeling et ahlbood ). In 
Figure [2]3 the existence of soft diffuse emission 
is emphasized in the context of the DSS op- 
tical image, where the diffuse X-ray emission 
nicely fills in the cavity of the HII region. This 
component will be discussed in Paper IV. 

2.1. Source Finding and Photon Event 
Extraction 

For identifying X-ray point sources, first we 
assemble a large number of candidate sources 
using a variety of techniques and criteria, in- 
cluding image reconstruction and visual in- 
spection. The source searching for NGC 2244 
is performed on the merged fields from Ob- 
sID 1874 and ObsID 3750. For each of the 
ACIS fields, twelve different images were cre- 
ated: soft (0.5-2 keV), hard (2-7 keV), and 
full (0.5-7 keV) X-ray wavebands with four 
different pixel binning scales (4x, ~2x, Ix, 
and 0.5 X a sky pixel). The wavdetect program 
(JFreeman et al.l l2002l ) was run with wavelet 
scales from 1 to 16, 8, 4, and 2 pixels in steps 
of \/2 (for the four different binnings respec- 
tively) and a source significance threshold of 
1 X 10^^ (which is very sensitive but permits 
some false sources) on each of the images de- 
scribed above. These twelve source lists were 
merged, with the source position from the 
highest-resolution image retained, to generate 
a single list of candidate sources. 

To take advantage of the sub-arcsecond 
point spread function (PSF) at positions 
around the aimpoi nt, we app l ied a subpixel 
positioning code (JMori et al.l I2OOII ) to im- 
prove spatial resolution in the inner part of 
the field. An image reconstruction with the 
Lucy- Richardson maximum likelihood algo- 
rithm ( Lucvlll974l ) was performed in the cen- 
tral 50" X 50" around HD 46150 (Figure E]) 
in ObsID 3750 (examples of maximum like- 
lih ood image reconstruc t ion can be foun d 



m 



Townslev et"ai] l2006al : IWang et al.l I2OO7I ). 



Eighteen additional candidate sources from 



the image reconstruction were added to the 
source Hst. Adaptive-kernel smoothed flux 
images in the three energy bands were also 
created with csmooth to help visually identify 
additional faint potential sources. Source lists 
from all ObsIDs (except the recent observa- 
tion, ObsID 8454) were then merged to form 
a master candidate detection list. 

The source finding procedure described 
above results in a total of 1452 potential 
sources identified for five ObsIDs (omitting 
the westernmost field ObsID 8454). A pre- 
liminary event extraction for the potential 
X-ray sources was made with our customized 
IDL scrip t A CIS Extract (v ersion 3.98; here- 
after ^K lBroos et al.ll2002l ). Using the AE- 
calculated probability Pb that the extracted 
events are solely due to Poisson fluctuations 
in the local background, source validity can 
be statistically evaluated while taking into 
account the large distorted PSFs at far ofF- 
axis locations and spatial variations in the 
background. After a careful review of the net 
counts distribution and IR counterparts fre- 
quency for all candidate sources, we rejected 
sources with Pb > 0.014, i.e. those with a 
1.4% or higher likelihood of being a back- 
ground fluctuation. The trimmed source list 
includes 1314 valid sources. 

Since our data analysis involves multiple 
Chandra observations, for convenience and 
to avoid repeating source designations we di- 
vide the X-ray point sources into NGC 2244 
sources and RMC sources based on positions. 
According to the stellar density distri butio n 
of 2MASS sources in NGC 2244 (OboO.'J'). 



the CO emi s sion m aps (|Williams et al.l 11995 
Hever et all l20nfil'). t he IRAS 60//m emis- 



sion ( Cox et al.l Il99d ) , and radio continuum 
(jCelnikI Il985l ). we assign all X-ray sources 
within 20 arcmin of the cluster central posi- 
tion (R.A.=06'*31™59.*9, Dec.=-F04°55'36") 



as potential NGC 2244 sources. This is indi- 
cated by the large circle in Figure [TJ The 
resulting cluster extent is consistent with 
the size of this massive open cluster deter- 
mined by systematic studies from the All- 
Sky Compiled Catalog of 2.5 Million Stars 



([Kharchenko et al.l |2005| ). We remind the 



reader that the dividing line is not unique; 
ambiguity of exact physical associations cer- 
tainly exists for sources located in the inter- 
face between the HII region and the molecular 
cloud to the east, and between the main NGC 
2244 cluster and secondary NGC 2237 cluster 
to the west. In this paper (Paper I) we focus 
on a total of 919 sources located within the 
NGC 2244 cluster region (as defined above); 
the rest are presented in Papers lu and III. 

The 919 valid sources are divided into a pri- 
mary list of 805 highly reliable sources (Pb < 
0.001; Table[2]) and a secondary list of 114 ten- 
tative sources with Pb > 0.001 likelihood of 
being spurious background fluctuations. Ta- 
bleland Table El have format s that are iden- 
tical to Tables 1 and 2 in iTownslev et al 



(|2006al l. IWang et al.l (|2007l 'l. and iBroos e~ 
20071 ). A detailed description of the table 



columns is given in the table footnotes. 

2.2. Source Variability 

One of the most notable characteristics 
of PMS stars is flaring in the X-ray band, 
and a few extraordinarily powerful X-ray 
flares hav e been reported in several PMS 
stars (e.g.. Ilmanishi et al.ll200ll: iGrossq et al 



20041: iFavataet al.ll2005l: ICetman et ahlbood : 



Wang et al.l l2007l : iBroos et al.l l2007|). A 
Kolmogorov-Smirnov (K-S) test is performed 
by AE on each observation in order to eval- 



^http : //www . astro . psu . edu/xray/docs/TARA/ | 



A small number of sources that are in ObsID 1875 and 
listed in Tables 2 and 3 here will be discussed scien- 
tifically in Paper II with the RMC sources. To avoid 
confusion, they will not be presented in RMC source 
list tables, but we will list these sources separately in 
Paper II for consistency between the papers. 



uate X-ray light curve variability within that 
observation, comparing the source event ar- 
rival times to that of a uniform light curve 
model. Seventy-eight sources display signifi- 
cant variability {Prs < 0.005 in column 15 of 
Tables [2] and [3]) and 14 of them have more 
than 200 net counts. Six of these highly vari- 
able light curves, for sources having more than 
500 counts, are shown in Figured 

Due to the short duration (20 ks) of most 
of our Rosette observations, few flares are ob- 
served in their entirety. However ObsID 3750 
is 75 ks long and completely captures a gi- 
ant flare from source 7^634 (average luminos- 
ity log Lt^c = 31.3 ergs s~^). The count rate 
during the flare peak is ~100 times that of 
the quiescent level. The estimated peak lu- 
minosity is log Lj^c = 32.3 ergs s~^. The 
shape of the flare is not symmetric, with a 
~1.5 hr rising phase and a ~4 hr decaying 
phase, resembling the fast rise and slow de- 
cay X-ray fl ares commonly seen in the COUP 
young stars ( Favata et al.ll2005l ). Its IR coun- 
terpart (§ [3]) does not show if-band excess, 
and the color and magnitude are consistent 
with a T Tauri star. The X-ray spectrum is 
very hard {kT ~ 7 keV). Multiple flares are 
seen in #691 with different intensities. Source 
#919 is only covered in the 20ks observation, 
but it shows a "flat-top" flare, characterized 
by a fast rise from low flux to high flux and a 
constant high flux level. 

2.3. Spectral Fitting 

For brighter sources with photometric sig- 
nificance Signify 2.0 (column 12 in Tables [2] 
and [21), the extracted spectra were fit us- 
ing single temperatur e and two-temper ature 
apec thermal plasmas (jSmith et al.ll200lh and 
power law models subjected to an absorbing 
column (A'^il/) of interstellar material with the 
XSPEu^ package (version 12.2. lap, lArnaud 



19961 ). based on source spectra, background 
spectra, ancillary response functions (ARFs) 
and redistribution matrix functions (RMFs) 
from ACIS Extract. The best-fit model was 
achieved by the maximum likelihood method 
( Cashlll979l ). Abundances of 0.3 Zq were as- 
sumed for the automated fitting performed by 
AE. 

In general, we prefer unconstrained fits (in- 
cluding power law fits) to constrained fits. 
The single temperature thermal plasma apec 
model is the default model used for spectral 
fitting. For sources brighter than 100 counts, 
if a one-temperature thermal plasma model 
did not fit the data well, a two-temperature 
thermal plasma model or variable abundance 
vapec thermal plasma model was invoked. A 
power law model was adopted if it represented 
the data more adequately than the thermal 
model (visually or with improved statistics) 
or if the thermal model required nonphysical 
parameters (e.g., kT S> 15 keV). Note that 
the adopted model should not be used to in- 
fer the nature of the source; a source best fit 
by a power law is not necessarily an AGN. 
For a number of sources that required a very 
hard thermal plasma and were identified with 
known stellar counterparts or exhibited flar- 
ing light curves that were suggestive of PMS 
stars, we truncated the plasma temperature at 
kT = 15 keV and adopt the thermal plasma 
model. When no model was acceptable, we 
froze the parameter kT = 2 keV in the ther- 
mal m odel, a typica l value for young PMS 



stars ( Getman et al.l l2005bl : iPreibisch et al. 



20051 ). and then fit for the absorbing column 



density Nh and the normalization parame- 
ters. The brightest source in the field is the O 
star HD 46150 (ACIS #373), which has 3588 
ACIS counts in 94 ks of exposure, or a count 
rate of 0.12 counts per CCD frame. This is 
not bright enough to cause multiple photon 
events in a single frame that would corrupt 



http://hea3arc.gsfc.nasa.gov/docs/software/lheasoft/xai IiHl0aglcf P^^*'^^^ ^^^i^g' ^^^^^ '^ ^oes not warrant 



corre ction for photon pile-up (JTownslev et al, 
20021 ). 



Spectral analysis results for the 630 sources 
with Signif > 2.0 are presented in Tables [H 
(588 sources; thermal plasma fits) and [5] (42 
sources; power-law fits). The table notes give 
detailed descriptions of the columns. Best-fit 
absorbing column densities range from neg- 
ligible to log Nh ~ 23.3 cm^^, equivalent 
to a visual a bsorp tion of Ay ~ 120 mag 



( Vuong et al.l 120031 ). Temperatures range 



from kT ~ 0.2 keV to the hardest sources 
truncated at kT = 15 keV. The range of ob- 
served total band (0.5 — 8 keV) absorption 
corrected luminosities able to be derived from 
spectral modeling is 29.1 < logL^^c ^ 32.4 
ergs s~^. Assuming a 2 keV plasma temper- 
ature and an average Ay = 1-5 mag visual 
extinction (logNfr ~ 21.4 cm~^ absorbing 
column), PIMMq^ gives an apparent total 
band luminosity logLj ~ 28.7 ergs s~^ es- 
timated for the faintest on-axis detection in 
Table El A conservative estimate for limiting 
sensitivity of the entire NGC 2244 observa- 
tion is log Li ~ 29.4 erg s~^ assuming an 
average extinction Ay = 1.5 and a 10 counts 
detection; the exact value depends on off-axis 
location and absorption. 

3. Identification of Stellar Counter- 
parts and Their Properties 

3.1. Stellar Counterparts Matching 
and IR diagrams 

We associate ACTS X-ray sources with op- 
tical and near IR (ONIR) sources using po- 
sitional coincide nce criteria, a s des cribed in 
the Appendix of iBroos et al.l toOJO § The 



''The Portable, Interactive Multi-Mission Simu- 
lator is software for high-energy astrophysicists, 
written and maintained by Koji Mukai. See 

[http: //heasarc.gsf c .nasa.gov/docs/software/tools/[ 

® Software implementing the matching algo- 
rithm is available in the TARA package at 



optical and infrared catalogs from recent lit- 
erature and observations that we adopted 
for counterparts identification include: UBV 
photometry of NGC 2244/Mon 0B2 (Massey 
et al. 1995 =MJD95), UBVIBa photome- 
try of NGC 2244 (Park & Sung 2002 =PS02), 
BVIRRa photometry of NGC 2244 (Bergh5fer 
& Christian 2002 =BC02), the Whole-Sky 
USNO-Bl.O Catalog (Monet et al. 2003 
=USNO), 2MASS All-Sky Catalog of Point 
Sources (Cutri et al. 2003 =2MASS), and the 
University of Florida FLAMINGOS Survey 
of Giant Molecular Cloudqj. The reference 
frame offsets between the ACIS fields (as- 
trometrically aligned to the Hipparcos frame 
using 2MASS sources in the data reduction) 
and the catalogs are 0.4" to MJD95, 0.3" to 
PS02, 0.3" to BC02, 0.2" to USNO, and 0.2" 
to FLAMINGOS. These offsets were apphed 
before matching sources. 

Likely associations between ACIS sources 
and ONIR sources are reported in Tabled 712 
of the 919 ACIS sources (77%) have an ONIR 
count erpart identified. Since iPark &: Sung 
(2002) provide the full list of their high pre- 
cision photometry, visual photometry is re- 
ported in the priority order of PS02, BC02, 
MJD95, and USNO. USNO photometry is 
photographic wit h ~ 0.3 magnitud e photo- 
metric accuracy (JMonet et al.l l2003l ) . JHK 
magnitudes from FLAMINGOS photometry 
are reported if available for Chandra sources. 
The FLAMINGOS photometric data was 



'http : //www . astro . psu . edu/xray/docs/TARA/ ' 
^PI: Elizabeth A. Lada; details about the 

FLAMINGOS Survey can be found at 

http : //flamingos . astro . uf 1 . edu/sfsurvey/sf su rvey .html [ 
The instrument design and performance of FLAMIN- 
GOS are described in Elston et al. (2003). An 
overview of the instrument is also available at 
http: //www.gemini . edu/sciops/instriiments/f la mingosl[ 
The FLAMINGOS observations of the Rosette Com- 
plex fields and IR data reduction are described in 
Roman-Ziifiiga (2006). Note that both K and K-short 
(Ks) filters are available in FLAMINGOS imaging 
mode, and the observations here are taken in K band. 



zero-pointed to 2MASS. (See Roman-Zuniga 
et al. 2007 for details on photometry. A 
photometric catalog for optical stars in NGC 
2244 is presented in RL07 and the complete 
catalog from the FLAMINGOS survey of the 
Rosette will be reported in Roman-Ziihiga 
et al. 2007.) For areas that were not cov- 
ered by the FLAMINGOS survey and for 
bright stars that were saturated (H < 11 
mag, Roman-Ziihiga 2006), 2MASS photom- 
etry is reported. The SIMBAD and VizieR 
catalog services are used for complementary 
information. Notes to other published charac- 
teristics of the selected sources can be found 
in the table footnotes. 

Figure shows the NIR J - H vs. H - K 
color-color diagram for 617 out of 919 Chan- 
dra stars with high-quality JHK photometry 
(error in both J — H and H — K colors < 0.1 
mag) listed in Table [H Most Chandra sources 
are located between the left two dashed lines, 
the color space associated with diskless young 
stars (Class III objects) that are reddened 
by interstellar extinction. We emphasize that 
while this region is usually filled with field star 
contaminants in NIR-only studies, here nearly 
all of these stars are cluster members (§3.2). A 
concentration of cluster members subjected to 
Ay ~ 1 — 2 mag (assuming late type stars) is 
apparent, centered at J—H=0.25, H—K=0.8. 
To the right of this reddening band are 38 
i^'-band excess sourceqj, defined as stars that 



^Note that the definition of _R'-band excess sUghtly 
varies when consider e d by different researchers. In 

iRoman-Zufiiga et al.l ((20071), a K-hand excess star 
is required to have colors that lie above J — H = 
0A7(H-K)+0A6 (locus of the Classical T Tauri Stars; 
CTTS), in addition to requiring that the star is located 
to the right of the reddening band for zero age main 
sequence dwarfs. The region below the CTTS locus 
could include detections of unresolved galaxies. But 
for the IR counterparts to our X-ray selected sample, 
this color space has little contamination from galaxies, 
and likely contains Herbig Ae/Be stars in the young 
cluster. Therefore no further constraint is applied in 
defining the if-excess. 



have colors (J-H) > 1.7{H-K)+2a{H-K). 
The presence of K-hand excess is frequently 
used to identify stars still possessing inner 
disks that are relatively hot (T ~1200K) and 
dusty. All except three of these stars occupy 
the color space between the middle and right- 
most dashed lines; they are likely PMS stars 
with circumstellar accretion disks (Class II 
objects). 

The three stars (#44, #564, and #805) lo- 
cated beyond the right-most dashed line show 
large color excess {E{H - K) = 0.87, 0.82, 
and 0.38, respectively). Stars in this domain 
are likely surrounded by extended envelopes, 
and hence a re classified as candidate Class 
I pro tostars ( Kenvon et al.lll993l : IStrom et al. 
19951 ). Sources #44 and #564 are particu- 



larly interesting {E{H — K) > 0.8). Source 
#44 was identified as a member of NGC 2244 
(jOgura k. Ishidalll98lh and assi gned spectral 
type B7Ve ( VerschuerenI Il99ll ). I t is fur 



ther c lassified as a Herbig Be star in iLi et al. 



(J2OO2I ) based on the tenuous nebula seen in 
their KPNO Hq, image, large [V-25/im] color, 
and its confirmed late B spectral type. Source 
#805 is very faint (J ~ 18.5, close to the 
FLAMINGOS imaging sensitivity limit) and 
may be a low mass protostar or an embed- 
ded background objec t with large H — K color 
(iFroebrich et al.ll2005l l. 

Figure [6] shows the NIR J vs. J — H color- 
magnitude diagram for the same stars shown 
in Figure [5j Known OB stars are located at 
the top and reddened from the ZAMS with 
Ay ^ 1 mag. Unhke NGC 6357 (jWang et al. 



20071) and other more heavily obscured clus- 



ters (jWolk et al.1 l2006l : iBroos et al.ll2007l l. we 
do not locate any new candidate massive stars 
earlier than BO. However two likely new BO- 
B2 stars (#900 and 919) and a dozen new late 
B candidates are found. 

Individual masses and reddening of the 
lower mass stars can be estimated assum- 
ing that they are reddened from the 2 Myr 



isochrone. Adopting a larger distance of d = 
1.6 kpc does not have a significant effect liere; 
the mass estimates for the low mass stars are 
only affected by ~0.1 Mq. The majority of 
ACIS stars appear to be concentrated around 
0.7 < J - iJ < 1 and 13 < J < 16, consistent 
with F, G, K, and M stars (0.1 < M < 2 M©) 
reddened by 1 < Ay < 2 mag. Those showing 
i^-band excess seem to have larger reddening 
compared to the Class III objects. Sources to 
the left of the 2 Myr isochrone are probably 
field stars older than the NGC 2244 popula- 
tion (see §3.2p . In addition, to the right of 
the ZAMS track, around twenty stars have 
inferred stellar masses < 0.1 M©. The IR 
photometry for these J > 17 mag stars may 
be less reliable, especially the fainter ones 
that are close to the detection limit. These 
are likely a mixture of the lowest mass clus- 
ter members, background stars, and a few 
IR luminous extragalactic sources. The high 
percentage of NIR excess so urces among the 



faintest stars may not be real; iFroebrich et al 



( 20051 ) have shown that background stars em- 
bedded in distant clouds have an overall larger 
H ~ K color. Contamination by faint back- 
ground IR sources is also discussed in Roman- 
Zuhiga (2006). 

The remaining ~200 ACIS sources with- 
out matched counterparts are likely to be a 
mixture of newly discovered embedded mem- 
bers of the Rosette complex and distant back- 
ground stars and extragalactic sources (see 
§3.21 for estimated fractions). Cumulative dis- 
tribution functions of the median photon en- 
ergy are shown in Figure [7| for sources that 
have matched counterparts and those that 
do not have counterparts. The harder me- 
dian photon energy of sources that do not 
have counterparts indicates that these sources 
are deeply embedded or behind the cloud 



{Ax 



> 



10 mag). The new low mass cluster 



stars reveal themselves because of strong X- 
ray flares due to magnetic activity. Such X- 



ray discovered stars are commonly found in 
the r nolecular clouds surrounding other clus- 
ters (ICetman et ahl lJoOSal: IWang et al.ll2007l : 
Broos et al.ll2007l ) . Some of these may be very 
young protostars with local absorptio n in an 
envelope or disk; ICetman et al.l ( 20071 ) found 
in IC 1396N that sources with logiV^ > 23.0 
cm~^ are protostars with dense envelopes. 
Their spatial distribution (Figure Ej) shows 
that they are distributed around the bright 
stars and along the rims of the Rosette neb- 
ula where infrared source identification may 
be challenging. 

3.2. Identification of Likely Contami- 
nants 

As emphasized in §1, stellar X-ray emission 
decays rapidly only after ~ 100 Myr and show 
only a small dep endence of X-ray luminosit;' 



[200; 



on disk accretio r i (IPr eibisch &: Feigelsoii 



Preibisch et al. 12005: Telleschi et al. 






20a , 
20m . 



Therefore X-ray surveys have very high effi- 
ciency in detecting disk-free PMS stars over 
MS stars, complementing traditional optical 
and infrared surveys of star forming regions. 
We evaluate the level of contamination by ex- 
tragalactic X-ray sources and Galactic disk 
stars following the simul ations descr i bed i n 
Getman et £01(120061 ^ and lWang etlH (|2007l ). 
The calculations indicate that < 35 sources 
are extragalactic (§ l3.2.ip . ~20 are foreground 
and ~16 are background field stars (§ I3.2.2p . 
These ~70 contaminants constitute ~ 8% of 
the 919 ACIS sources. 

3.2.1. Extragalactic Contaminants 

Following Getman et al. (2006), we 
construct Monte Carlo simulations of the 
extragalactic population by placing arti- 
ficial sources randomly across the detec- 
tor. We draw incident fiuxes from the X- 
ray background log A^ — log S distribution 



(JMoretti et al.l l2003l ). and power law pho- 



ton indices for the sources are assigned con- 
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sistent with flux-de pendencies described by 
Brandt et alj (J200lh . The spectrum of each 
simulated source was passed through a uni- 
form absorbing column density log A'^^^ ~ 21.9 
cm~'^, the HI column density through the en- 
tire G alactic disk in the direc tion of NGC 
2244 (|Dickev fc LockmanI Il990l ^. After ap- 
plying local background levels found in our 
ACIS image, we calculate the photometric 
significance of each fake source and then re- 
ject the weak extragalactic sources that would 
have fallen below our source detection thresh- 
old. The simulations suggest that r^ 80 ex- 
tragalactic sources may be detected in our 
ACIS-I field and ~ 35 may have photometric 
significance Signif > 2.0 (column 12 in Ta- 
ble [2|). The true number is probably smaller 
as we did not account for the patchy distri- 
bution of molecular cloud material. The best 
candidates for the extragalactic contaminants 
are those sources whose X-ray spectra are 
best fit by a power law, do not have bright 
ONIR counterparts, and do not display the 
characteristic PMS X-ray flares. 

3.2.2. Galactic Stellar Contamination 

Monte Carlo simulations of the Galactic 
stellar population expected in the direction of 
our ACIS fields (/ = 206.3, b = -2.1) were ex- 
amined, based on the stellar population syn- 
thesis model of Robin et al. (2003; henceforth 
the Besancon modeljj. In addition to the 
smooth absorption component provided in the 
model, we added a local absorption at d = 1.4 
kpc with a low Ay ~ 1 mag, as inferred from 
the average NIR reddening of Chandra sources 
shown in Figure [6j Within the ACIS field-of- 
view, the Besangon model predicts ~ 850 fore- 
ground MS stars {d < 1400 pc) and ~ 4300 
background MS stars, giants, and sub-giants 



within the FLAMINGOS imaging sensitivity 
limit around J < 19 mag. 

About 25% of the foreground stars in the 
Besangon simulation have ages less than 1 
Gyr, the younger population that could be de- 
tected in X-ray surveys. Most reside at dis- 
tances from 0.8 to 1.4 kpc; ~54% are M stars, 
~28% are K stars, ~13% are G stars, and 
3% are more massive. Following Getman et 
al. (2006), we convolved the Besangon model 
populations with the X-ray luminosity func- 
tions of stars in the solar neigh borhood mea- 



sured from ROSAT survevs (ISchmitt et al. 



1995; Schmitt 1997: Hiinsch et al. 



199 



i). Lu- 



minosities were adjusted to account for the 
different ROSAT and Chandra spectral bands 
follo wing the stellar har dness-luminosity rela- 
tion (JGiidel et al.lll998l ). Following the same 
procedure we used for extragalactic sources 
( ^3.2.ip . we applied our ACIS detection pro- 
cess to these simulated field stars. A typical 
Monte Carlo run predicts that 16 — 20 fore- 
ground stars will be detected in the ACIS-I 
exposures of NGC 2244. 

Five optical counterparts to our X-ray 
stars are identified as foreground stars based 
on their known spectral types and positions 
in photometric diagrams. These are #321, 
429, 583, 627, and 67 8 . Proper-rnotion data 
(JMarschall et al.l Il982l : bias et al.l l200fil ) also 
suggest that they have low probability of be- 
ing member stars of the cluster. In the NIR 
color magnitude diagram, there are ~ 30 
stars located between 2 Myr isochrone and 
the ZAMS track. Those close to the 2 
Myr isochrones are still consistent with be- 
ing members due to uncertainty in the age 
and distance of the cluster. About fifteen 
stars close to the MS track are more likely to 
be unrelated foregrou nd main sequence stars 



( Getman et al.l 120061 ). Available X-ray ab- 



* These calcul ations are made wit h the Web service pro- 
vided by the R obin et al.l (|2003j group at Besangon at 
[http: //bison, obs-besancon.fr/modele 



sorption columns derived from spectral fits 
mostly require log Nh = 20.0 cm~^, sugges- 
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tive of being foreground candidatecl. Their 
probabilities of being members derived from 
proper motion data range between 0% and 
80%, which are less conclusive. With the 
above five confirmed foreground stars, alto- 
gether twenty optical starq^ are noted as the 
best candidates for foreground stars in the Ta- 
ble [6] footnotes. This number is quite consis- 
tent with the Besancon simulated foreground 
population. 

The Besangon model is again convolved 
with X-ray luminosity functions (XLFs) to 
simulate the number of stars behind the 
Rosette star forming region that may enter 
our X-ray sample. The model predicts ~ 18% 
F dwarfs, ~ 41% G dwarfs, ~ 33% K dwarfs, 
~ 1%M dwarfs, and ~ 6% giants. We use the 
dwarf XL Fs established for the solar neigh - 
borhood (jSchmitt et al.l Il995l : ISchmittI Il997l ) 
and adopt the XLF for gia n ts obt ained from 



Table 2 of Pizzolato et al. 



(|200d l. Typical 
runs of simulations result in ~11 dwarfs and 
~5 giants that are detectable in our Chandra 
observation. To the left of the ZAMS track in 
the CMD and for J > 18, there are 16 stars^l] 
whose locations match the Besangon model 
predicted background population of MS stars, 
subgiants, and giants, which we note as best 
candidates for background stars in Table O 



^Spectral fits to low counts sources can give Nh val- 
ues with large uncertainties, hence we do not use low 
Nh to identify foreground sources. Low absorption 
columns are considered supporting evidence for fore- 
ground stars suggested by their locations in the NIR 
color magnitude diagram. The spatial distributions of 
the low Nh sources with high counts (> 30 counts) 
and with low counts (< 30 counts) are similarly dis- 
persed through the field, consistent with being part of 
the foreground population. 

^° They are associated with Chandra sources #49, 67, 
96, 105, 129, 163, 321, 429, 552, 583, 603, 627, 678, 
743, 757, 785, 839, 854, 859, and 889. 

"They are #68, 157, 158, 181, 226, 253, 323, 325, 347, 
409, 442, 583, 609, 747, 835, and 880. 



4. Global Properties of the Stellar 
Cluster 

4.1. X-ray Luminosity Function and 
Initial Mass Function 



As noted in iFeigelson et al.l (J2005l ). the 



XLF (which is directly measured here) can 
be considered to be the convolution of the 
IMF (which is unknown) and the X-ray-Mass 
Luminosity (L^ — M) correlation (which is 
measured in the COUP studies; Preibisch 
et al. 2005). Using the best-studied Orion 
Nebula Cluster XLF (COUP XLF) and IMF 
as a calibrator, the NGC 2244 XLF can be 
used to probe the IMF of the stellar clus- 
ter and to estimate the total X-ray emitting 
population. Such a population analysis has 
been made for Cep 0B3b, NGC 6357 and 
M17 dCetman et al.l I2OO6I : IWang et al.l l2007l : 
Broos et aDl2007l ). In the following XLF anal- 
ysis, we use the hard band XLF rather than 
the total band XLFs, since the unknown soft 
component of heavily absorbed X-ray sources 
can introduce a large uncertainty in both the 
observed and absorption corrected total band 
X-ray luminosity. 

By counting the number of sources in dif- 
ferent X-ray luminosity bins, we construct 
the absorption corrected hard band (2-8 keV) 
XLF {Lh^c) for all unobscured NGC 2244 X- 
ray sources {MedE < 2.0 keV) with derived 
X-ray luminosities in Figure [9^. Here we ex- 
clude the five known foreground stars and OB 
stars with spectral types earlier than B3 to be 
consistent with the Orion cool stars sample. 
The absorption corrected hard band fluxes 
{Fh^c) derived from XSPEC spectral fitting 
are used to obtain luminosities assuming a dis- 
tance of 1.4 kpc. As the template, we also 
show the XLF of the COUP unobscured pop- 
ulation (839 cool stars; Feigelson et al. 2005). 

The NGC 2244 XLF is largely consistent 
with the COUP XLF, suggesting an NGC 
2244 population comparable to the ONC. 
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However, at log L/i^c ^ 29.8 ergs s ^ the 
slope of the NGC 2244 XLF seems steeper 
than the COUP XLF. At logL^^^ = 30.9 ergs 
s~^, the NGC 2244 bin is short 9 stars com- 
pared to the ONC XLF; the luminosity bins 



at log Lhc = 29.5 - 29.7 ergs s" 



are 



40 



stars higher than the ONC XLF. Admittedly 
these deviations from the ONC XLF are not 
significant {2a), more closely following the 
shape of the ONC than the Cep B or the 
M17 XLF does (see Figure [9b). The apparent 
steeper slope is not an artifact of our distance 
estimate or detection completeness limit. We 
comfortably detect 15 count sources at any 
off-axis location; the X-ray luminosity of the 
detected sources is roughly logL^ ~ 29.4 ergs 
s^^ (corresponding to ^^ 0.5Mq from the 
Lx — M relation) and L^^c ~ 29.2 ergs s~^. 
The luminosity bins with the steep slope seen 
in the NGC 2244 XLF are much brighter than 
the completeness limit. 

A few possible extrinsic reasons may ac- 
count for the slope deviation from the ONC 
XLF. First, while the log-normal COUP 
XLF represents the best data and provides 
a good observational template, the underly- 
ing physics and the variations from it among 
clusters remain to be explored. There is ev- 
idence that the XLF may n ot be identical 



in al l regions (see review by iFeigelson et al. 
20071 ). Figure Eb suggests that the XLFs of 



Cep B and M17 may vary from the ONC XLF 
in different ways. The L^ — M correlation is 
derived from the very young ONC, but may 
vary for the older clusters like NGC 2244. 
Second, the exclusion of some hot stars to 
construct an XLF comparable to the COUP 
XLF could implicitly affect the resulting XLF. 
We examined luminosities (logL/^c) of these 
OB stars, which mostly contribute to the lu- 
minosity bins around logL^c ~ 30.0 ergs s~^; 
an XLF including them still shows the ob- 
served deviation from the ONC. 

There are two possibilities intrinsic to the 



NGC 2244 cluster that may be responsible for 
the apparent steeper slope: (i) The NGC 2244 
population is the same as the COUP popula- 
tion, but there is an excess of ~ 50 stars in 
the luminosity range 30.0 < logL/j,, ^ 30.4 
(solar mass stars as inferred from the Lx — M 
relation); or (ii) The NGC 2244 population 
is ~ 1.2 times larger than the COUP pop- 
ulation, but NGC 2244 is ~ 20 stars defi- 
cient in stars with 30.4 < logL/^^c ^ 31.0 (in- 
termediate mass stars). In either case, the 
NGC 2244 XLF deviates from a scaled ver- 
sion of the COUP XLF in a manner similar 
to that seen in th e Cep B/0B3b field studied 
bv lGetmanet al.l (|200fil ). 

To improve the statistics, remove possi- 
ble binning effects, and further investigate 
how the NGC 2244 XLF compares to that 
of the ONC and other clusters, we derive 
the cumulative distribution of X-ray lumi- 
nosities for the unobscured NGC 2244 X-ray 
sources as well as the unob scured COUP, 
Cep B/0B3b (age~l-3 Mvr. ICetman et al 



2006 ). NGC 6357 (age~l Myr, 



2007), and M17 (age~l Myr, 



Wang et al 



Broos et al. 



20071 ) populations. The resulting cumula- 



tive XLFs are shown in Figure [9b. To avoid 
confusion caused by incompleteness in de- 
tections, we cut off the cumulative XLFs 
at the corresponding completeness limit for 
each region (see Figure [9]). At a given X- 
ray luminosity above the completeness limit, 
the ratio between the cumulative numbers 
of sources from two populations reflects the 
relative scaling between the two populations 
that are more luminous than this limit. As 
a result, the unobscured population of NGC 
6357, M17, and Cep B, is ~ 5 times, ~ 3 
times, and ~ 0.4 times of the size of unob- 
scured ONC population, respectively. These 
are consistent with previously reported val- 
ueso The cumulative X-ray luminosity func- 



In NGC 6357, log Lh,c is ~0.5 dex higher than the 
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tion of NGC 2244 closely follows the COUP 
XLF, although the deviation found in Figure 8 
(deficit at 30.4 < log Lh^c ^ 31.0 and excess at 
log-L/jc ^ 30.4) can still be clearly seen. De- 
pending on the treatment of this deviation, 
the NGC 2244 population ranges from 1.0 to 
1.2 times the ONC population. 

Further investigations were made to exam- 
ine the possible excess of sources with 30.0 < 
log Lh^c ^ 30.4 ergs s^^. To identify a previ- 
ously unknown cluster in the field that may 
contribute extra stars to NGC 2244, we in- 
spected the spatial distribution of sources that 
have luminosities in the excess bins, but did 
not find any apparent clustering. To test 
whether this excess comes from contamina- 
tion of X-ray bright non-members, we re- 
moved candidate contaminants (foreground 
stars, background stars, extragalactic sources) 
as suggested in § 13.21 and reconstructed the 
XLF. The number drop mainly appears in low 
luminosity bins while the excess is still signifi- 
cant at 30.0 < Lh,c ^ 30.4 ergs s~^. Therefore 
we conclude that excluding candidate contam- 
inants would not alter the XLF, since they 
do not contribute much to the high luminos- 
ity bins that characterize the X-ray emitting 
population. 

4.2. Initial Mass Function and K-band 
Luminosity Function 

To examine whether the deviation in the 
XLF is a reflection of an intrinsic anomalous 
IMF in the NGC 2244 cluster, we perform 
two tests on the IMFs using NIR data. One 
experiment is to use the location of X-ray 
stars in the IR color magnitude diagram to 



other clusters because candidate OB stars are included 
in the sample (|Wang: et al.ll2007l ). Note that, although 
the unobscured population of NGC 6357 is larger than 
that of M17, the obscured population in Ml 7 is sig- 
nifica ntly larger than that of NGC 6357 (jBroos et al.l 
120071 ). which makes the estimated total populations of 
the two clusters comparable. 



derive their masses and construct an approx- 
imate IMF. The exact mass will not be as ac- 
curate as measured from spectral types and 
an HR diagram, but their statistical distribu- 
tion should be sufficient for our interest here. 
If there is indeed an excess of X-ray sources 
with logLhc ~ 30.0 — 30.4 ergs s~^ (option i 
i n j^4.1), from the ern pirical L^ — M relation 
(jPreibisch et al.ll2005l ). we would expect to see 
an excess of stars around a solar mass. Fig- 
ure [TO] shows the IMF constructed from NIR 
estimated masses for unobscured COUP stars 
and for NGC 2244 stars. No excess is appar- 
ent for stars in the solar niass ra nge, after the 
ONC IMF (JMuench et al.l I2OO2I I is scaled to 
match the NGC 2244 IMF. Instead, a deficit 
of intermediate mass stars around 2 — 3 Mq in 
the X-ray selected sample is apparent (option 
ii in §4.1). 

A second test is to obtain a statistical 
sample of cluster meinbers and construct 
a KLF (jLada fc Ladal Il995l ) following the 



Cep B /0B3b KLF analysis in iGetman et al. 
(|2006l l. We used 2MASS Ks data to con- 



stuct the KLF since it was spatially com- 
plete toward NGC 2244 and it is photometri- 
cally complete to a similar mass limit as our 
Chandra data for the Rosette. A control field 
from 2MASS centered at {a, 5) =(6''30™00^ 
-h3°45'00") [J2000] is used for background 
population subtraction (extinction toward 
this control field and in the foreground of 
NGC 2244 is low; see Li 2005 for its 2MASS 
color-color diagram). The resulting KLF is 
similar to that derived by Li (2005) with a 
power law slope {dlog N {K g) / dK s) ~ 0.3. To 
get the IMF, we convert Kg magnitude to 
logM using the 2 Myr theoretical isochrone. 
The relation between K^ and log M frorn 



the theoretical isochrone ( Siess et al.l l200d ) 



can be approximated well with a power law 



as demonstrated by ICetman et al.l ( 20061 ) 



The resulting IMF derived from i^-band 
star counting (dot-dashed line in Figure [TO]) 
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is consistent with the IMF estimated from 
NIR properties of Chandra sources (solid his- 
togram in Figure [To]) . Note that the mass de- 
rived using K magnitude as a proxy for pho- 
tospheric luminosity can be biased towards 
higher mass for K-excess sources. However, 
our IMF estimated from KLF is justified; the 
fraction of iT-excess sources in our sample 
is low and the IMF inferred from the KLF 
is coarsely binned (Alog(M/MQ) ~ 0.2) to 
account for the effect of K-excess on the esti- 
mated mass. 

This KLF suggests that our X-ray sample 
is largely complete down to ~ 0.5Mq, con- 
sistent with the mass limit estimated from 
the X-ray completeness limit. It also sug- 
gests that our X-ray selected sample is miss- 
ing ~ 20 intermediate-mass stars around 3 
Mq. It is not surprising that the detection 
efficiency of intermediate-mass stars in the X- 
ray band is not as deep as in optical or IR 
(jSchmitt et al.lll985l ). The X-ray production 
mechanism is not well understood for stars in 
the intermediate- mass range (^ I5.2p . In many 
cases of X-ray detections of intermediate-mass 
stars, the X-ray emission in fact comes from 
a low mass companion to the intermediate 
mass star. Compared to the ONC, some 
of the intermediate-mass stars in NGC 2244 
may be X-ray quiet because of the absence 
of low- mass companions. This is reflected as 
the deficit of log Lh,c ~ 30.4 — 31.0 ergs s~^ 
sources in the XLF. The dynamical evolu- 
tion history of the young cluster might result 
in a lower fraction of intermediate-mass stars 
with secondaries. Surveys of intermediate- 
mass member stars in NGC 2244 to mea- 
sure their multiplicity, such as the binarity 
study of Sco 0B2 in t erme diate-mass stars by 
Kouwenhoven et al.l (J2007l ). are needed. Un- 
der this assumption, the "excess" of log Lh,c ~ 
30.0 — 30.4 ergs s~^ sources would then be the 
result of a slightly larger population of NGC 
2244 stars compared to that of the ONC. 



Based on all the above analysis, we con- 
clude that the unobscured X-ray emitting 
NGC 2244 population is about 1.2 times 
larger than the known unobscured population 
in ONC, or ~1000 stars with logLt > 27.0 
ergs s~^ (the COUP detection limit). The 
obscured population, estimated from sim- 
ilar XLF scaling, is ~ 500 stars. Given 
that the total COUP sample accounts for 
~ 75% of the ONIR sarnple of the ONC in 
Hillenbrand &: HartmannI ( 19981 ) . we estimate 



that the total size of the NGC 2244 stellar 
population is around 2000. |lj (|2005l ) gives 
a census of ~1900 NGC 2244 members es- 
timated from the spatially complete 2MASS 
analysis. 

Although the massive end of the stellar 
complement in NGC 2244 was well-known 
from early studies (e.g.. lOgura &: Ishidalll98ll ). 
the l ow-mass pop ulations were poorly identi- 
fied. IPerea ( 199ll ) suggested that there may 
not be M < 4 Mq stars in the NGC 2244 
cluster. Perez (1991), Massey et al. (1995), 
and Park & Sung (2002) investigated the IMF 
for NGC 2244 in optical, and reported a top 
heavy IMF with a fiat power law slope F 
(dlog7V(logM)/(ilogM) ~ -0.7, although 
Park & Sung (2002) are cautious due to the 
incompleteness of their intermediate- and low- 
mass population. Indeed, their optical sam- 
ple becomes incomplete for stars with masses 
lower than ~ 3Mq. In comparison to the 
traditional IMF studies, our high sensitivity 
X-ray sample of stars largely benefits from our 
robust membership criteria and our ability to 
identify the low-mass members in a reliable 
manner. The IMF slope obtained from the 
X-ray selected membership, which is nearly 
complete to O.5M0, gives F ~ —1.1 instead of 
F ~ —0.7, consistent with the Orion IMF. 
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4.3. Spatial Structure of the Stellar 
Cluster 

The morphology of young clusters, includ- 
ing dependency on stellar mass, provides 
clues for cluster formation and dynamical 
evolution. As morphological studies based 
on optical or infrared samples are compli- 
cated by patchy extinction, nebular contam- 
ination, confusion with field stars, and bias 
towards stars retaining protoplanetary disks, 
the spatial distributions of X-ray identified 
stars in populous young clusters should be 
excellent laboratories to explore their ori- 
gins and dynamical evolution. For example, 
an aspherical shape or clumpy distribution 
would reflect unequilibrated initial conditions 
while a spherical shape with mass segregation 
would indicate we ll-developed virialization 
(jClarke et al.l boool ;). In the Orion A cloud, 
the ONC, NGC 2024, and associated molec- 
ular filaments have flatte ned shapes (jLada 



199ll : iFeigelson et all l2005l ) which have been 



attributed to global gravitational collapse 
of an elongated cloud (JHartmann &: Burkert 
20071 ). In contrast, the rich NGC 6357 and 
M 17 clusters appear spherical but with sub- 
clusters that may reflect distin ct (perhaps 



triggered) subcluster f ormation (IWang et al 



20071 : iBroos et al.ll2007l ) . The absence of mass 
segregation can either reflect a young stellar 
system that has not yet achieved dynami- 
cal relaxation, or a mature system where 
many of its massive members have been 
ejected by few-body interactions in the core 
( Pfiamm- Altenburg fc Kroupall2006l ). 

It was recognized in the 2MASS study by Li 
(2005) that the apparent center of the large- 
scale annulus (see Figure lb) defining the op- 
tical Rosette Nebula at (a, 6) = (6 ^31'"56^ 



04°59'56") (jOeura k Ishidal Il98lh is offset 
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from the center of the IR surface density dis- 
tribution at (a, 6) = (6''31™59.9^ -F04°55'36"). 
He interprets this offset as a projection effect 
where the Rosette Nebula resembles a tilted 



cylindrical cavity in th e molecular cloud (see 
models in ICelnikI 119861 ). The distribution of 
the Class II sources in a rec ent Spitzer /IRAC - 
MIPS survey of NGC 2244 (JBalog et al.ll2007l ) 
shows good agreement with the 2MASS data. 
Here our Chandra data show that the stellar 
concentration around HD 46150 is offset from 
the larger stellar distribution, irrespective of 
its relation to interstellar matter. The off- 
center massive star HD 46223 will be further 
discussed in § 14. 4[ 

A 'center' must be defined for structural 
analysis. In NGC 2244, the highest concen- 
tration of X-ray stars appears around the sec- 
ond most massive star, HD 46150 (05V), at 
(a, (5) = (6'*31™55^+04°56'34"). We there- 
fore treat HD 46150 to be the center for our 
radial profile analysis of this high density re- 
gion in § 14.3.21 although we are aware of the 
asymmetrical distribution of X-ray stars as 
noted above: HD 46150 lies about 2' north- 
west of the IR center defined in the 2MASS 
study (Li 2005). 

4-3.1. Morphology and Substructures 

Figure [11] shows a smoothed map of the 
stellar surface density for 572 lightly-obscured 
(median photon energy < 2keV) NGC 2244 X- 
ray stars. This map of smoothed spatial dis- 



(|2007l l. 
has been 



tribution is constructed f o llowin g IWang et al. 
hooii ) and iBroos et al 



A similar 
applied to 



smoothing technique 
2MASS data to identify large-scale struc- 
ture of the entire Rosette Complex (JLj l2005l ; 
Li fc Smithir2005al lbll. A ~ 20' x 20' grid is cre- 
ated to cover the stellar positions, and at each 
position the total number of sources within a 
0.5 arcmin radius sampling kernel is counted 
to estimate the smoothed stellar density. Only 
sources covered by both ObsID 1874 and Ob- 
sID 3750 are considered to guarantee roughly 
equal X-ray sensitivity throughout the field. 
The heavily-obscured sources are omitted be- 
cause many of them are expected to be back- 
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ground AGNs. 

The cluster shows an approximately spheri- 
cal structure that extends 8' (3.2 pc) in diame- 
ter, centered at (a, 5)=(06'^31™59^04°55'30")• 
This center, as well as the large scale structure 
and substructure seen in the X-ray-sampled 
cluster (Figure [TT|) . are in good agreement 
with the results derived from the surface den- 
sity of K-excess stars in the FLA MINGOS 



study ( Roman-Zuhiga et al.1 120071 ). It also 



matches the c enter defined from the 2MASS 
star-count (JLi| 120051 ) and the ce nter defined 



frorn the Spitzer Class II sources (jBalog et al. 
20071 ). The large-scale north-south asymme- 



try can be attributed to the off-center place- 
ment of HD 46150. The primary concentra- 
tion is seen around HD 46150. Five of these 
stars were noted by Sharpless (1954) as a 
visual compact subcluster, but we find ~50 
stars extending to a radius of 1' around this 
massive star. The central stellar surface den- 
sity here is ~700 stars per pc^; recall that this 
value is restricted to stars with masses above 
~ 0.5 M0 due to X-ray sensitivity limits. 

A secondary density enhancement of ~ 15 
X-ray sources (about a 3a enhancement) is 
seen 3' south of HD 46150 at {a, (5)=(6''31™56^ 
Different sampling scales are tested and this 
substructure persists for smoothing kernels 
with radii < 1'. The local density peak here 
has six stars tightly clustered within 20"; it is 
also apparent in the optical Ha ima ge and th e 
2MASS-KS image (Figure [H] andQ l2005l ). 
Assuming they are lightly-obscured late-type 
stars, the reddening indicated by their NIR 
colors is Ay ~ 1 mag, similar to that of 
NGC 2244 cluster. Their NIR estimated spec- 
tral types range from F to M if they are lo- 
cated at the same distance as NGC 2244. 

The existence of both substructures, around 
HD 46150 and 3' to the south, is direct ev- 
idence that the NGC 2244 cluster has not 
attained dynamical equilibrium. But perhaps 
most remarkable is the absence of companions 



around the most massive cluster member, HD 
46223 (04V). Nine X-ray sources lie within 
1' of HD 46223 compared to ~50 around HD 
46150. One possible explanation for the iso- 
lation of HD 46223 is that it was ejected by 
dynamical interactions within the HD 46150 
subcluster. However, it does not exhibit high 
proper motion (§4.4) and it seems unlikely 
that such a massive member, rather than less 
massive members, would be ejected at high 
velocity. We note, however, that an 04 su- 
pergiant has been reported to be probably 
ejected from Cyg 0B2 (Comeron &: Pasquali 
2007). Radial velocity measurement of HD 
46223 will be valuable to evaluate the ejec- 
tion scenario. 

The spatial distribution of the NGC 2244 
sources exhibits some common characteris- 
tics and notable di fferences when co mpared 
to M17 (see §3.1 in iBroos et al.ll2007l l. They 
both show the highest density of stars close 
to massive O stars and the concentration is 
largely spherical. However, the concentration 
in M17 is around the known early O stars in 
NGC 6618, while in NGC 2244 the concentra- 
tion is not around its massive star with the 
-|-04°54^4^*^?t spectral type (HD 46223) but with an- 
other O star, HD 46150. Both clusters show 
distinct substructures: in M17 an obscured 
small cluster is found (M17-X), which is seen 
as an elongation of the central cluster in M17; 
an unobscured substructure is also found in 
NGC 2244. M17 shows a triggered stellar pop- 
ulation along the shock front and the eastern 
edge of the M17-SW molecular core (south 
bar in Jiang et al. 2002). To the southeast 
of NGC 2244, a sequence of triggered embed- 
ded clusters also exists along the midplane of 
the RMC (Phelps & Lada 1997). The differ- 
ent appearances of the triggered populations 
are related to the geometric configuration of 
the dense giant molecular clouds relative to 
the HII regions. The "V"-shaped M17 is 
edge-on blister HII region emerging from the 
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surrounding molecular materials, while NGC 
2244 is located in an expanding HII bubble at 
the tip of the elongated RMC. 

4-3.2. Radial Density Profile 

The radial density profile for the NGC 2244 
cluster, centered at the stellar density peak 
around HD 46150, is shown in Figure [T3] with 
comparison profiles fro m optical/NIR and X- 
ray s t udies of the ONC (IHillenbrand &: Hartmann 
19981 : iFeigelson et all l2005l ) and from our 
X-ray study of P ismis 24 in NGC 6357 



19981: lAdams et alJ l200l|). ISchilbach et al 



(|Wang et al J 120071 ') . The radial profile of 



NGC 2244 has two distinctive components: 
a powerlaw structure around HD 46150 ex- 
tending 1.5', and a structure with a flat core 
and steeper dropoff extending from 1.5' to 
8'. The power law structure is centered on, 
but is much more extended than, the ~ 20" 
X-ray resolved subcluster around HD 46150 
shown in Figure 3. NGC 6357, and perhaps 
the ONC, have a similar radial profile with 
approximately the same powerlaw slope (Fig- 
ure USD. 

Figure O shows the inferred radial den- 
sity profiles in physical size units (par sees) for 
the three clusters, where the star densities are 
scaled to their estimated true densities based 
on the comparison of the XLFs shown in Fig- 
ure 9b. The stellar density of NGC 2244 has 
been scaled to 1.2 times the ONC population 
(§ ED, and NGC 6357 to 5 times the ONC 
population (IWang et al.l 120071 ). Omitting the 
central r~^ powerlaw structures, the profiles 
of these two clusters can be fit as isothermal 
spheres with core radii Tc = 1.2pc and r^ = 1.4 
pc, respectively. 

4.4. Mass Segregation 

The concentration of massive cluster mem- 
bers at the center and lower mass members 
at larger radii from the cluster center is com- 

mon ly observed in rich young star clusters 

(e.g. lCarpenter et al.lll997l : lHillenbrand &: Hartmann 



( 20061 ) investigate mass segregation in over 
600 open clusters with a wide range of ages. 
For their youngest clusters with ages ~5 Myr, 
some show mass segregation whereas others 
do not. 

Mass segregation can occur as a nat- 
ural consequence of dynamical relaxation. 
Details of the process have been debated. 
For the ONC Trapezium, some researchers 
argue that the dense collection of ~ 10 
OB stars is an imprint o f initial conditions 
(IBinnev &: Tremaind 119871 : iBonnell &: Davies 
19981 ). while others argue that the core has 
collapse d and many OB stars have bee n 
ejected ( Pflamm- Altenburg &: Kroupal l2006l ) . 



McMillan et al.l (|2007l ) suggest a model of se- 
quenti al mergers of mass se gregated subclus 



ters. 



Bonatto et aD (|2006l ) found that the 
M 16 cluster, with age ~1.3 Myr, has an over- 
all relaxation timescale around ~ 20 Myr, yet 
shows some degree of mass segregation at this 
young age. 

In NGC 2244 in the Rosette Nebula, the O 
stars are not highly concentrated, as shown in 
Figure [T5b (see also Figure la). The earliest 
0-type (04V) star in the cluster, HD 46223, 
has a rather puzzling location in the cluster 
near the southeast boundary of the nebula. 
Its proper motion is fairly sn iall {(in = —0 .2 



mas yr~ , /i^ = 0.4 mas vr~ : IZacharias et al. 
20041 ) which does not suggest ejection from 
the cluster center. The issue of mass segrega- 
tion has not been extensively investigated for 
this 2 Myr old cluster in the literature, mainly 
because the low-mass population was not ad- 
equately identified. With the proper motions 
and membership probabilities of stars in the 
NGC 2244 region derived fr om photographic 
plate data, IChen et al.l (J2007l ) suggest that the 
cluster shows evidence of mass segregation. 

Figure [T5b shows the cumulative radial dis- 
tributions for the massive stars with NIR- 
estimated masses M > 8M0 and for the X- 
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ray identified low mass stars (M < 2Mq). 
The distributions appear very similar, and 
a Kolmogorov-Smirnov test of the two dis- 
tributions does not show significant differ- 
ence. Thus, mass segregation is not present 
in NGC 2244 within our FOV. 

We estimate a two-body dynamical re- 
laxation time trpjax for the NGC 2244 clus- 



ter (e.g. 



Bonatto et Zl [ioO^) 



i^relax ~ 



{N/8lnN) X tcross where tcross 
is the characteristic crossing time for a star 
to travel through the cluster with radius 
R and velocity dispersion Vdisp- Adopting 
i? ~ 4 pc from the full cluster extent in 
Figure [TH a rough estimate for the unmea- 
sured velocity dispersion Vdixp ~ 3 km s~^ 
teinnev fc Tremainj [iQStI '). and N ~ 1900 
stars (§ 14.11) . we obtain treiax ~ 30 Myr for 
NGC 2244. As the cluster age is < 10% of 
this relaxation time, no significant mass seg- 
regation is expected from two-body dynam- 
ical interactions. If we consider only stars 
within the estimated core radius re = 1.2 pc, 
then treiax ~ 9 Myr which is still considerably 
larger than the age of the cluster. 

The absence of mass segregation is thus 
consistent with standard dynamical theory, 
and implies that NGC 2244 (unlike some other 
clusters) was not formed with a central con- 
centration of massive stars. The main chal- 
lenge for explaining the dynamical state of 
NGC 2244 is the difference between the dom- 
inant member HD 46150, which has a rich 
compact subcluster, and HD 46223, which is 
mostly isolated. 



Similar to NGC 2244, iBroos et al.l (J2007l ) 
show that many other massive stars are scat- 
tered all over the ACIS-I field in M17, in addi- 
tion to the concentration in NGC 6618. ACIS 
source ^51 is one of the most massive stars in 
the field, yet it sits > 6'from the center of the 
cluster. For comparison, HD 46223 is about 
~ 7' south of HD 46150 and ~ 5'south of the 
cluster center in NGC 2244. Not all massive 



stars are participating in the mass segregation 
seen in NGC 6618. One possible explanation 
is that the O stars are not all co-eval. Indeed, 
several mas sive protostars have been found in 



M17 (e.g., iNielbock et al.l l200ll : IChini et al 



20041 . I2OO5I ). Together with the presence of 



an UCHII region there, these young massive 
stars establish that the massive populations 
are not all the same age in Ml 7. Massive 
stars to the east might be older and belong to 
a wider OB association, while source #51, for 
example, might be younger. This could also 
be the case in NGC 2244. The late-0 stars 
are scattered as in M17. HD 46223 may be 
younger and not part of the same population 
as NGC 2244's central cluster. 

4.5. X-ray Stars with Infrared Excess 
Disks 

X-ray selected samples have several advan- 
tages over optical and IR samples (see review 
by Feigelson et al. 2007). X-ray emission 
arises from stellar magnetic activity which is 
enhanced 10^ — 10^ above main-sequence levels 
for stars during the entire age range of interest 
(< 0.1 to > 10 Myr), thus X-ray surveys suffer 
only a small number of field and extragalactic 
contaminants (~ 8% in observations of rich 
massive clusters), which are usually identifi- 
able (§3.2). They naturally deliver a nearly 
disk-unbiased sample of young starq^. The 
main disadvantage of typical X-ray surveys is 
their incompleteness in detecting the lowest 



^■^There may be additional complications: X-ray selected 
PMS samples suffer a small bias against accreting stars 
in the 0.5 — 8 keV band because Class II systems are 
on average ~ 2 times fainter than Class III systems 
(Preibisch et al. 2005; Telleschi et al. 2007), and there 
may also be a small bias toward accretion systems in 
the soft < 1 keV band due to emission at the accre- 
tion shock. However accretion variation s do not cause 
X-ray variations in the Chandra band IjStassun et al.1 
|200y) . These are minor effects considering that the X- 
ray luminosity function spans 28 < log Lx < 32 ergs 
s^^ and it is dominated by flare emission. See discus- 
sion in Feigelson et al. (2007). 
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mass objects, which can be identified in high- 
sensitivity IR images. The complementary 
nature of the Chandra and IR data will pro- 
vide the best census to date for the young stel- 
lar population of NGC 2244. It is worth not- 
ing that because the different disk indicators 
in the IR trace thermal emission from circum- 
stellar materials of different temperatures, the 
inferred disk fraction of a young cluster gen- 
erally increases w hen measured in longer IR 
wavele ngths (e.g.. iLada et al.ll2004l ). For ex- 
ample, iHaisch et al.l (|2001al ) show that JHK 
observations alone are not sensitive enough to 
detect circumstellar disks in a complete and 
unambiguous manner. They find that 21% 
of the sources in the IC 348 cluster have K- 
band excess disks; the disk fraction rises to 
65% when using JHKL IR-excess emission. 
In this subsection we focus on the X-ray se- 
lected sample of stars with i^-band excesses. 

The spatial distribution of 38 X-ray stars 
with K-h&nd excesses attributed to inner pro- 
toplanetary disks (§ 13. 1|) is shown in Fig- 
ure dni Ten of these young stars cluster 
around HD 46150 (05V), and a few are 
around HD 259135 (B0.5V). 

A deficit of color-excess stars in the north- 
ern part of the nebula is seen: none of them 
are located in the northern part of the clus- 
ter where HD 46149 (08.5V) is located. This 
region is also devoid of dust e mission in the 
IRAS image (ICox et all Il990l ). The deficit 
of disky stars cannot be attributed to pho- 
toevaporation by OB stellar UV radiation 
and winds since a grouping of them is found 
around HD 46150. We consider two expla. - 
nations for this asymmetry. First, |lJ ( 20051 ) 
suggested that the Rosette Nebula is open 30° 
north of the line of sight. As gas and dust 
stream away from th e HII region (as in M 17, 
Townslev et al.ll2003l ). the stars in this region 
may have undergone a faster inner disk dissi- 
pation so that their disks no longer show K- 
band excess. Second, the star formation in 



NGC 2244 may have proceeded over a consid- 
erable time span along the north-south direc- 
tion, with the older population in the north- 
ern region. Similar spatial-age patterns have 
been f ound in other young clus ters (e.g., Cep 
0B3b, iBurningham et al.ll2005l '). 

For the ~2 Myr old NGC 2244 cluster, us- 
ing our X-ray selected sample with high JHK 
photometric quality (largely complete to 0.5 
M0), we derive an overall disk frequency of 
~6% for Chandra stars with mass M > O.5M0 
(assuming a presence of ~20 foreground field 
stars). The disk fraction is ~10% for stars 
with mass M > 2.0Mq and ~5% for stars 
with mass O.5M0 < M < 2.OM0. Using a 
large Chandra sample si milar to that studied 
here, IWang et al.l (J2007l ) reported a low frac- 
tion of /fg-band excess among intermediate- 
mass stars in the young massive star form- 
ing region NGC 6357 (~1 Myr old) and a 
similar result is reported for M 17 (~1 Myr 
old) in iBroos et al.l (J2007l ). These are consis- 
tent with the findings that optically thick cir- 
cumstellar disks are already rare among the 
intermediate-mass PMS stars with ages less 
than a few Myr and suggest that the disks 
are short lived for the massive stars ( e.g., 
Hillenbrand et al.lll993l : lNatta et al.ll200d l. It 



has been suggested that the disks around ear- 
lier type stars may evolve faster than around 
later type stars based on studies of the IR- 
excess fraction as a function of spectral type in 



a few clusters (ILada et alj 
2001bl : iLada et al.ll2006l l. 



M 



2OOOI: Haisch et al 



O (J2OO5I I used the 2MASS NIR sample 



after background population subtraction to 
derive a X^-excess disk fraction of ~ 20.5% 
above mass ~ O.8M0 for NGC 2244. The dis- 
crepancy between our results reflects the dif- 
ferent criteria in selecting excess sources: their 
K^-excess sources included 2MASS sources 
with large photometric error. With high qual- 
ity FLAMINGOS JHK photometry data, 
Roman-Zuhiga et al.l ( 20071 ) derive a lower K- 
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excess fraction of ~ 10% for K < 15.75 stars, 
although it remains shghtly higher than the 
IR-excess fraction in our X-ray-selected sam- 
ple. Their sample covers a larger mass range, 
probably d own to 0.1-0.2 Mm depend ing on 
extinction (|Roman-Zuhiga et al.l 120071 ). 

Balog et al.(2007) present a Spitzer survey 
of NGC 2244 covering 3.6 fun to 24 fim and 
estimate that the overall disk fraction in the 
cluster is 44.5% using the IRAC data (and 
MIPS data if available). The total cluster 
population in their field of view (0.5 square de- 
gree), including the cluster members without 
disks, is estimated by subtracting the aver- 
age background population, whereas our sam- 
ple identifies the diskless stars from their el- 
evated X-ray emission. But the discrepancy 
between our results is mainly due to the fact 
that Spitzer mid-IR observations are much 
more sensitive to circumstellar dust than K- 
band excesses. The overall IR-excess frac- 
tion among NGC 2244 Chandra-selected stars 
will be studied with a deep Spitzer survey 
(PI: Bouwman) and compared to the IR- 
only determined excess fraction, aiming to 
evaluate the different samples of young stars 
selected through IR colors and those iden- 
tified in X-rays. A similar study has al- 
ready taken place in the Serpens cloud core 
([Winston et al.l 120071 ) , which analyzed com- 
bined Spitzer and Chandra observations of the 
embedded stellar cluster. 

The overall low X-excess disk frequency 
seen here (~ 6% for Chandra stars with mass 
M > O.5M0; ~ 10% for FLAMINGOS stars 
with mass M > O.IMq), at a cluster age of 
2 Myr, may imply a faster disk dissipation 
while the cluster is immersed in the hostile 
environment of UV radiation and strong stel- 
lar winds of many massive stars. In clusters 
of similar age but without the presence of 
O stars, the K-band excess fractions appear 
higher. F or example, in the 2 .3 Myr old IC 
348 cluster iLada fc Ladal (llQQSl ^ find that 20% 



of the stars (M > O.OSMq, estimated from 
the X-magnitude completeness limit) have K- 
excess disks. The recent Spitzer observations 
of IC 348 find a total disk frequency of ~ 50% 
among stars with mass M > O.IMq based 
on infrared excess between 3.6 and 8 .0 fiui 
200fil : iMuench et aP l2007l ^. not 



(Lada et al. 



significantly higher than the ~ 45 % disk frac- 



tion i n NGC 2244 reported in iBalog et al 



( 20071 ). This implies that although the inner 
disks around the NGC 2244 stars may dissi- 
pate faster, there seems no difference in the 
overall disk dissipation. 

However, iBalog et al.l (J2007l ) did notice a 
lower disk fraction (27%) among stars close 
to the NGC 2244 O stars (separation d < 
0.5 pc), hinting a faster disk evolution near 
the massive stars. A number of observational 
and theoretical studies have already demon- 
strated the photoevap oration of disks by ex- 



ternal radiation (e.g. , lO'Dell fc Wone 
Johnstone et al. 1998: Hollenbach et al. 



1996 



200d; 



Throop fc Ballvl |2005|). IBalog et al.l (120061 ) 
presented 24/im images of three protoplane- 
tary disks being photoevaporated around high 
mass O stars, including one disk close to the 
05 star HD 46150 in NGC 2244 with an es- 
timated mass loss rate lO"^'^ — 1O~^M0 yr~^. 
We do not detect X-ray emission from the IR 
point source in this cometary structure, likely 
due to its low mass. Based on the Spitzer iden- 
tified Class II and Class I sources, it is further 
suggested that the effect of massive stars on 
the circumstellar disks is significa nt in the im- 



medi ate vicinity of the hot sta rs ([Balog et al. 



20071). JGuarcello et al.l (|2007|) reported evi 



dence that the spatial distribution of the stars 
with a circumstellar disk in NGC 6611 is anti- 
correlated with the distribution of OB stars. 
Their findings also suggest that UV radiation 
from OB stars does have an impact on the 
evolution of the disks close to massive stars. 
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5. X-rays across the Mass Spectrum 
5.1. X-rays from Massive Stars 

One of the important discoveries of early 
Einstein observations was tlie soft X-ray 
emission from individ ual early type O stars 
( Harnden et al.l Il979l ). Most 0-type stars 



were found to be soft X-ray emitters {kT < 1 
keV) with X-ray luminosities L^ ~ 10^^ — 10^'^ 
erg s"^. A canonical relation between X- 
ray luminosity and bolometric luminosity 
of Lx/Lhoi ~ 10~^ was proposed and con- 
firme d from Einstein and ROSAT observa- 
tions ( Pallavicini et al.lll981 : Chlebowski et al. 



198^). lBergh6fer et alJ (jl997|) extended the 
same relation down to stars of later spec- 
tral type (B1~B1.5). Recent Chandra stud- 
ies of 07-B3 stars in Orion (COUP) found 
both a soft wind-emission component and a 
hard flaring component in many OB stars, 
and a larger dispersion was found for late O 
and early B stars (log(L,r/-L/)n;) ~ —4 to —8; 
Stelzer et al.ll2005l ). However, when only con- 
sider iiig_^£2lM_§,niission in the 0.5-2.5 keV 
band, ISana et al.l ( 20061 ) derived a tight scal- 
ing law log (L^/Lboi) = — 6.91±0.15 for 0-type 
stars with a deep XMM-Newton observation 
of NGC 6231. 

Wind-shock models were developed to ex- 
plain the X-ray emission from massive stars, 
where small-scale instabilities in radiatively- 
driven stella r winds from massive stars pro- 
duce s hocks (Lucv &: Whitelll98d : IOwocki et al. 



19881 : lOwocki fc CohenI Il999l l To account 



for the observed X-ray emission line pro- 
files and hard, variable continuum eniission 
fe.g..lCorcoran et al1ll994l:lEvans et al.lf2004l : 



the hard components in Chandra gratings 
spectra of 6^ Ori C (06V). In some cases, 
the anomalously hard and luminous X-ray 
component {kT > 10 keV and Lh ~ 10'^^ 
ergs s~^) implies close binarity, as power- 
ful winds in two massive components collide 
and shock to produce very h i gh energy X- 
rays (e.g.. Pollock et al.l 12005 : Skinner et al.l 



20061 : iBroos et al.ll2007l ). ISchulz et al.l (120061 ) 
observed a large X-ray outburst in 9'^ Ori 
A, which can be attributed to reconnection 
events from magnetic interactions between 
the binary stars. 

Due to its richness in population, NGC 
2244 offers an excellent opportunity to study 
X-ray emission in OB stars. Table [7] sum- 
marizes the detection/non-detection of O and 
early B-type stars in our observation, along 
with their optical, IR, and X-ray properties. 
Chandra spectra are shown in Figure [l7l We 
detected all 9 OB stars with spectral types 
BO. 5 or earlier that were in the field. The 
two early O stars in NGC 2244 exhibit soft 
{kT < 1 keV) and strong {L^ ~ 10^^ erg 
s~^) X-ray emission as expected in the clas- 
sical wind-microshock regime. All of the 2.3 
Myr old NGC 2244 O stars show soft X-ray 
emission {kT < 1 keV), which is similar to 
the soft thermal spectra seen in most 0-type 
stars in other star forming regions (e.g., Rauw 
et al. 2002, Rho et al. 2004, Skinner et al. 
2005, Sana et al. 2006, Albacete-Colombo et 
al. 2007). Nevertheless unusual cases of O 
stars characterized by a hard spectrum have 
been reported in a number of observations 
(e.g., Gagne et al. 2005, Broos et al. 2007, 
Tsujimoto et al. 2007). Based on analysis of 
X-ray grating spectra of the Orion Trapezium 
stars, Schulz et al. (2003) proposed that the 
presence of hard X-ray emission indicative of 



Waldron et al.ll2004l : IStelzer et al.ll200,5l ). more 

complex models were invoked such as the 

magne tically channeled wirid shoc k (MCWS) 

model teabel fc Montmeri^|l997al lbl: fud-Doula fc OwJ^^Slil coronal plasma {kT > 1 keV) indicates 

the presence of coronal magnetic fields. Lin- 
sky et al. (2007) compared the X-ray proper- 
ties of young massive stars with known mag- 



2OO2I ). iGagne et al.l (|2005l ) shows that the 



MCWS model with strong line-driven winds 
can adequately reproduce both the soft and 
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netic fields (e.g., theta^ Ori C, 06V, t ~ 0.3 
Myr; M16 ESI) with older stars with weak or 
no magnetic fields (r CMa, 09.51, t < 12 Myr; 
^ Ori, 091, i ~ 3 — 5 Myr). The very young 
stars with magnetic fields show high coronal 
plasma temperatures {kT ~ 2 — 5 keV), which 
may be heated as in the MOWS model (Gagne 
et al. 2005), and the total X-ray luminosi- 
ties are far larger than the expected X-ray 
luminosities from micro-shocks in the stellar 
winds. The absence of such a hot compo- 
nent in the X-ray spectra of the older massive 
stars (including our NGC 2244 O stars) might 
imply a timescale of the presence of coronal 
magnetic fields in the massive stars; the hot 
component may be restricted to stars much 
younger than 2 Myrs. A detailed analysis of 
X-ray properties of a large sample of O stars 
will be required to link the magnetic fields in 
massive stars to the cluster ages. 

We detected 6 out of 14 B stars with spec- 
tral types B1~B3. This low X-ray detection 
rate among B stars is consistent with other re- 
cent Chandra and XMM-Newton observations 



of ni a ssive star for i ning r egions (IWang et al.l 



2007; Broos et al. 2007; Sana et al. 



200i 



m- 



The X-ray emission from early B stars is con- 
sistently harder than that from the O stars 
(higher kT in Table [7]), which suggests that 
unseen late-type companion star s rather than 



the B star itself is responsible (jStelzer et al, 
200i). 



The earliest exciting star in this com- 
plex, HP 46223, is of spectral type 04V 



(Walborn et al. 



2OO2I ). Its X-ray spectrum 
is adequately fit by a soft kT = 0.3 keV sin- 
gle temperature plasma subjected to Nh = 



4 X 



10^^ cm- 



absorption (Figure [TTjl . HD 
46150 (05V) is the visually brightest early 
type star in the cluster. A two temper- 
ature plasma model fit {kTi = 0.2 keV, 
kT2 = 0.6 keV) is needed to describe the 
X-ray spectrum, with an absorption column 
Nh = 2.5 X 10^^ cm-2. Their X-ray lumi- 



nosities are similar, Lt^c ~ 2.5 x 10^^ ergs 
s~^. The light curves of the 0-stars are ex- 
amined and no variability is suggested by the 
K-S statistics. In contrast, both 04V stars 
in Ml 7 show a very hard plasma compo- 
nent [kT > 10 keV) and are nearly an order 
of magnitude brighter in their intrinsic full- 
band X-ray emission ( Broos et al.l 120071 ). If 
this hard emission is caused by close bina- 
rity or fossil magn etic fields, as suggested by 
Broos et al.l ( 20071 ). then it is likely that the 



early O stars in NGC 2244 lack at least one 
of these features; their soft X-ray emission 
suggests that these are single O stars without 
close, massive companions and/or that they 
do not possess strong magnetic fields. 



Berghofer &: ChristianI (J2002l ) reported the 



X-ray and optical luminosities of NGC 2244 
early type stars from ROSAT PSPC and 
HRI observations and BVI photometry, and 
concluded that they are consistent with the 
canonical relation Lx/Li^oi ~ 10~^. The 
Lx/Lboi relation for OB stars in NGC 2244 de- 
termined from our Chandra data is shown in 
Figure [TBI Statistical tests for the correlation 
between Lx and L^oi were perform ed using the 



ASU RV survival analysis package (llsobe et al 



19861 ). To take into account the available up 



per limits, the generalized Kendall's tau corre- 
lation test for censored data is adopted. The 
null hypothesis (a correlation is not present) 
probability is P < 0.01%, supporting a sig- 
nificant Lx — Lfjoi correlation. As shown in 
Figure M as well as Table H the NGC 2244 
O stars closely follow the X-ray to bolometric 
luminosities Lx/Lboi ~ 10~^ ratio, although 
the B spectral type stars show larger (yet still 
< 0.5 dex) scatter. For comparison, previ- 
ously reported Lx vs. L^oi values for addi- 
tional OB stars from the massive star forrn - 
ing regions Orion (1 Myr; Stelzer et al.ll2005l ). 



NGC 6357 (1 Mvr: IWang et al. 



200 



and 



M17 (1 Mvr jBroos et al.ll2007l )are also shown 
in Figure [m The overall scatter in Lx/Lf,oi is 
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considerably smaller for O stars from different 
clusters (~2 orders of magnitude) than for B 
stars (~4 orders of magnitude). 

5.2. X-rays from Intermediate Mass 
Stars 



X-ray emission from intermediate-mass 
stars with spectral types mid-B to A is un- 
expected since no X-ray production mecha- 
nism is known; they lac k strong stellar winds 
and convective surfaces (Berghofer Sz Schmit 



ral variability {Pks ^ 0.005), which strongly 
supports emission from flaring, possibly from 
unresolved low mass companions. For exam- 
ple, source #804 shows a big flare during our 
long observation (Figure H]). The rising and 
decay times during the flare are rather sym- 
metric. The count rate during the peak of 
the flare is 18 times higher than that of the 
quiescent level. As noted in § 13.11 source #44 
is a previo usly identified Herbig Be star in 



NGC 2244 (iLiet aLlliooJ l. We detected 40 



1994 ; iBerghofer et al.lll997l : iStelzer et al.ll200: 

20061 ) . However X-ray detections of PMS 

intermediate-mass stars known as Herbig 



net counts at its optical position in 94 ks. The 
fit to its X-ray spectrum indicates a low ab- 
sorption and a hard {kT = 2.8 keV) plasma 



Ae/B e stars are widely repor t ed fe.g..lZinnecker &: Preibiscl 




iDamiani et al 



Stelzer et al. 



a] nonvariable lightcurve. 
n *= 



1993: ^Berghofer et al 



2005 : Hamaguchi et al. 



i). A systemat ic Chandra archival study 
( Stelzer et al.ll2006l ) rules out radiative winds 
as the origin for X-ray emission in HAeBes 
based on the observed high X-ray tempera- 
tures. Thus an X-ray generating mechanism 
from magnetic flares similar to late type stars 
or emission from an unknown/unresolved 
companion is favored, although the role of ac- 
cretion f or the production of X - rays remains 



uncle ar (iHamaguchi et al.l l2005l : IStelzer et al, 
20061 ). 



Using masses estimated from the NIR color 
magnitude diagram (Figure 6), we detect 
around 50 stars in the intermediate mass 
range 2Mq < M < 8Mq. The absorption- 
corrected X-ray luminosities of the detected 
sources in the 0.5-8.0 keV band are in the 
range of logLj^c ~ 29.5 — 31.8 ergs s~^. 
This is fully consistent with the level of 
X-ray emission detected in HAeBes from 



Hamaguchi et al.l (|2005l ) and IStelzer et al 



(J2006l ). but also with the level expected from 
late-type companions. Optical spectral classi- 
fication of these stars may further clarify the 
link between the X-ray emi ssion propertie s 



and the spectral types (e.g., iLi et al.l |2002| ) 



5.3. X-rays from Other Interesting 
Sources 

Herbig-Haro Jets and Knots Two 

optical jet systems, namely Rosette HH 1 
(R.A.=06'^32"^20.'*76, Dec=04°53'02.9") and 
HH 2 {R.A.=06^ 32"" U.'U, Dec.=05°02'17.95" 
[J2000]), ha .ve been discovered in th e R osette 
Nebula (HboOSl: iLi k Rector] bood: Hl2005l: 



Meaburn et al.l l2005l : iLi et al.l I2OO6I I. HH 1 



Eight of these stars show significant tempo- 



consists of a collimated jet originating from 
a faint optical star with a mass-loss rate 
M ~ 1O~^M0. We have detected X-ray emis- 
sion (ACIS source #743) coincident with the 
location of the exciting source of HH 1, a 
weak-lined T Tauri star (F8Ve, Li et al. 2005). 
At the location of the base, knot, or termi- 
nal shock of the jet, the distribution of X-ray 
photons is consistent with background counts. 
No concentration of photons is coincident with 
the shock structure seen in optical near the 
end of the collimated jet. The X-ray spec- 
trum of the star is soft, with kT ~ 0.8 keV 
and a negligible absorption column. It has 
been suggested that the combination of low 
extinction and high Lyman photon flux inside 
the Ro sette Nebula rnakes the jet optically 
visible (JMeaburn et al.ll2005l ). 

A group of bright ionized knots in the 
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Rosette Nebula were proposed to be colli- 
sionally ionized, either by bow shocks formed 
around globules by the strong winds from O 
stars which are then overrun by an expand- 
ing shell, or collim ated flows of shocked ga s 
driven by the wind ( Meaburn Sz Walshlll986l ). 
Chen et al.1 (|2004l ) found an X-ray source in 
the nebulous region D with high speed knots 
and identified two stars as counterparts within 
ROSAT positional errors. We examined our 
Chandra image with much higher resolution 
at all knot locations. No X-ray emission was 
found to be associated with these high-speed 
knots, but we detected point sources coinci- 
dent with optical stars embedded in the neb- 
ulous region. ACIS source #678 (= Chen et 
al. source 30) matches the position of HD 
259210, a likely foreground star with spectral 
type AlV. 

Binaries— The eclipsing binary V578 Mon, 
used to determine the age and distance to 
NGC 2244, is detected in our X-ray image 
with 186 net counts (^^476). It consists of two 
early B-type stars, one of t he very few massive 
eclip sing systems known (JHarries &: Hilditch 
19981 ). Its orbital period is precisely deter- 
mined in the optical to be P = 2 .40848 ± 
0.00001 days (JHensberge et al.ll200d ). Using 
the eclipse ephemeris and the date at the be- 
ginning of our observation, we examined the 
X-ray light curve together with the optical 
light curve with orbital phases (Figure [20]) . 
No significant variability is suggested by a 
K-S test for the X-ray light curve, although 
dips might be seen in the X-ray light curve 
around phase 4> = 0.85 and (f) = 1.0, where 
the primary eclipse is expected. This could 
be one of the rare cases where X-ray eclipses 
can be used to constra. i n ern itting geometry 
(e.g.. ISchmitt &: Favatalll999l ). However this 
could simply be a statistical fluctuation given 
the limited number of counts. The spectral fit 
gives a plasma temperature of kT ~ 1.6 keV 
with low absorption, logNn = 20.6 cm~^. 



Park &: Sund ( 20021 ) noted a suspected 
PMS binary system ([PS02] 125 and [PS02] 126). 
An optical spectrum of the unresolved bi- 
nary shows Hn in emission a nd Lil 6708A in 



absorption ( Chen et al.ll2004l ). confirming its 
youth. We resolve the system as a close pair of 
ACIS sources (#243 k #242). Their spectral 
fits give a rather hard kT ~ 2.9 keV for #243 
and kT ^ 1.5 keV for #242; they share the 
same logA'^^ = 21.4 cm~^. The light curve of 
#242 is constant while #243 shows possible 
but not significant variation {Pks = 0.03). 

Magnetic Star- ISagnulo et~all (|2004l ) 
discovered an extraordinarily strong mag- 
netic field in the very young cluster member 
NGC 2244-334 (=[0181] 334; spectral type 
B3; R.A.=06'^32™51.^79, Dec=-h04°47'l6.l" 
[J2000]), ranking as the second strongest lon- 
gitudinal field known among non-degenerate 
stars (after HP 21 5 441, or Babcock's star; 
Borra fc Landstreetl Il978l ). We detected 
a 10 count X-ray source (#899, CXOU 
J063251.79-F044715.9) at its position in the 
20 ks observation (it was not covered by the 
deep observation). 

Elephant Trunks— ISchneps et al.l (1980) 
identified several spectacular elephant trunk 
globules in the northwest part of the Rosette 
Nebula. Only one of them, a small isolated 
globule denoted Rl, is in our field of view. 
One X-ray source (#169), probably by chance 
superposition, is located 6 arcsec away from 
the bright rim of this dark globule. No other 
X-ray/IR source can be found inside the glob- 
ule. 

In the southeast quadrant of the nebula 
towards the RMC, another molecular pil- 
lar is prominent. This region i s also high 



lighte d in a recent Spitzer survey (JBalog et al 



20071 ). as the size of the pillar is compara- 



ble to the l a rgest "pillar-of-creation" in M16. 



Chen et al.l (|2004l ') noted shocked gas near the 



pillar, perhaps due to strong winds from a 
star nearby that was matched to one of the 
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ROSAT sources. In the vicinity of this pil- 
lar, a luminous X-ray source (#919, CXOU 
J063309.61-F044624.3) is detected in our ob- 
servation, but it is not located at the tip of 
the elephant trunk. Its X-ray spectrum can be 
fit well with log Nh = 21.2 cm~-^ and a hard 
(fcT = 3.3 keV) plasma. With an unusually 
high luminosity of log Lf^c ~ 32.2 ergs s~^, 
it is comparable to the earliest O stars. The 
light curve is variable, as shown in Figure [H 
The count rate doubles after the first 6 ks 
and remains in a high state for ~13 ks. Its 
IR counterpart is also bright, with a i^-band 
magnitude of 9.6 mag that is similar to the 
observed B0-B2 stars in the field. Its location 
in the color-magnitude diagram also suggests 
a spectral type of BO-Bl. No ET-excess is 
seen. As an early B-type star, its X-ray vari- 
ability can be explained by an unresolved late 
type companion, although the lightcurve does 
not follow the typical PMS fast-rise and slow- 
decay phase. A similar transition between 
high and low states seen in 9 ^ Ori A is investi- 
gated by ISchulz et al.l (120061 ) and interpreted 
as possibly the reconnection from magnetic 
interactions in a close binary system. 

6. Summary 

We present a high spatial resolution X-ray 
study of the NGC 2244 cluster in the Rosette 
Nebula obtained via deep Chandra observa- 
tions. Our main findings follow: 

1. We detect 919 X-ray sources with a lim- 
iting X-ray sensitivity of Lt^c ~ 1 x 10^^ ergs 
s~^. Positional coincidence matching yields 
a total of 712 ONIR counterparts. We esti- 
mate 8% extragalactic and Galactic contami- 
nation. The rest of the X-ray sources without 
ONIR counterparts are hkely new NGC 2244 
members clustered around the massive star 
HD 46150 or deeply embedded in the cloud. 
The X-ray detected population provides the 
first deep probe of the rich low mass popula- 
tion in this massive cluster. 



2. The locations of most ACIS sources 
in the color-magnitude plot indicate a large 
population of 2 Myr old PMS low mass stars 
(M < 2M0) subject to a visual extinction of 
1 < Ay < 2 at 1.4 kpc. We derive an overall 
-fC-excess disk frequency of ~6% for stars with 
mass M > 0.5Mq using the X-ray selected 
sample, slightly lower than the 10% if -excess 
disk fraction using a FLAMINGOS selected 
sample. Both fractions are significantly lower 
than the 45% mid-IR disk fraction in a Spitzer 
sample that is more sensitive to disks. We em- 
phasize that the combination of young stars 
identified in X-rays (mostly Class III stars) 
and those selected through IR colors will pro- 
vide the best census to date for the young stel- 
lar population of this region. Three objects 
have Class I colors. 

3. The derived XLF {Lh,c) for NGC 2244 
is compared to the XLFs of the ONC, M17, 
Cep B, and NGC 6357; this exercise indicates 
that the unobscured population in NGC 2244 
is 1.2 times larger than that of the ONC, or 
~ 1000 stars detectable in COUP-sensitivity 
X-ray observations. Taking into account the 
obscured population, the total stellar popula- 
tion in NGC 2244 is ~ 2000, in good agree- 
ment with the estimated population from the 
spatially complete 2MASS study. The XLF 
and KLF suggest a normal Salpeter IMF for 
NGC 2244; we do not confirm a top heavy 
IMF reported from earlier optical studies. 

4. We examine the spatial distribution of 
the X-ray identified NGC 2244 cluster mem- 
bers; the stellar surface density map suggests 
a spherical cluster with substructure. We con- 
firm the existence of a sub cluster around HD 
46150 with ~ 50 members in a 1 pc region; a 
second small subcluster consisting of a num- 
ber of late type stars is also found. The 
04 star HD 46223 has few companions. The 
radial density profile of NGC 2244 shows a 
larger relaxed structure around the central 
subcluster. Similar structure is seen in NGC 
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6357. No evidence for significant mass seg- 
regation is found in this cluster. Altogether 
we suggest that this 2 Myr cluster is not dy- 
namically evolved and has a complex star for- 
mation history. Our results will strongly con- 
strain models of the cluster formation process. 

5. We detected all 9 OB stars with spec- 
tral types BO. 5 or earlier, but only 6 out of 
14 B stars with spectral types B1-B3 in our 
field of view. X-ray spectra for the massive 
stars in NGC 2244 all show soft emission. We 
confirm the long-standing log(L^/Lf,o;) ~ — 7 
relation for the NGC 2244 O stars. Large scat- 
ter around this correlation was found for the 
B stars. 

6. We report X-ray emission detected from 
a few interesting individual objects, including 
the ionizing source of the optical jet Rosette 
HHl, binary systems, a magnetic star, and a 
possible X-ray luminous uncataloged massive 
star. 
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Table 1 
Log of Chandra OBSERVATIONS 



Target 






Obs ID 


Start Time 
(UT) 


Exp 


osure Time 


Ain 


ipoint 


Roll Angle 
(deg) 


Modc^ 




a.I2000 


'^J2000 




Rosette Field 1... 
Rosette Field 2... 
Rosette Field 3 






1874 
1875 
1876 
1877 
3750 
8454 


2001 Jan 05 11:53 
2001 Jan 05 17:46 
2001 Jan 05 23:28 
2001 Jan 06 05:10 
2004 Jan 01 02:20 
2007 Feb 09 02:25 




19700 
19500 
19410 
19510 
75000 
20480 


06 31 52.85 
06 32 40.84 
06 33 17.15 
06 34 17.34 
06 31 56.45 
06 30 50.40 


+04 55 42.0 
+04 42 45.0 
+04 34 42.0 
+04 27 45.9 
+04 56 25.4 
+04 59 34.0 


335.86 
335.90 
335.76 
335.85 
351.87 
286.00 


F 
F 
F 


Rosette Field 4 

Rosette Nebula/NGC 2244 
NGC 2244 Satellite Cluster 


F 
VF 
VF 



^Thc observing mode: F— Faint, VF— Very Faint. 

Note. — Units of right ascension are hours, minutes, and seconds; units of declination are degrees, arcminutcs, and arcscconds. Exposure 
times are the net usable times after various filtering steps are applied in the data reduction process. The aimpoints and roll angles are 
obtained from the satellite aspect solution before astrometric correction is applied. 
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Table 2 
Chandra Main Catalog: Basic Source Properties 







Source 




Position 








Extr 


acted Counts 








Charact 


eristics 








ycq 


CXOU J 


".J2000 


^^.72000 


lirr 





Net 


ANct 


Bkgd 


Net 


PSt' 


Bignit 


logics 


Anom 


Var 


Ett'Exp 


^median 




# 




(dcg) 


(dcg) 


(") 


(') 


Full 


Full 


Full 


Hard 


Frac 










(ks) 


(keV) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


(17) 




1 


063114.36+045303.0 


97.809835 


4.884189 


0.9 


9.9 


23.3 


5.7 


3.7 


5.8 


0.91 


3.7 


<-5 


g''' 




15.0 


1.4 




2 


063117.04+045228.3 


97.821033 


4.874535 


0.7 


9.5 


32.0 


6.5 


4.0 


1.4 


0.90 


4.5 


<-5 




a 


17.0 


1.4 




3 


063118.11+045511.8 


97.825464 


4.919953 


1.0 


8.7 


12.4 


4.4 


2.6 


3.5 


0.90 


2.5 


<-5 




a 


16.3 


1.7 




4 


063118.29+045223.1 


97.826245 


4.873106 


1.0 


9.2 


14.6 


4.8 


3.4 


0.8 


0.89 


2.7 


<-5 




a 


16.8 


1.2 




5 


063118.76+045207.6 


97.828205 


4.868779 


1.1 


9.2 


9.7 


4.1 


3.3 


1.9 


0.89 


2.1 


-4.4 




b 


17.0 


1.1 




8 


063120.89+045003.8 


97.837054 


4.834396 


0.2 


10.7 


599.7 


25.8 


33.3 


32.1 


0.87 


22.8 


<-5 




a 


81.8 


1.0 




311 


063152.54+050159.1 


97.968929 


5.033106 


0.1 


5.8 


299.3 


17.8 


0.7 


0.6 


0.43 


16.3 


<-5 




a 


85.7 


0.9 




373 


063155.51+045634.2 


97.981330 


4.942835 


0.0 


0.4 


3588.5 


60.4 


0.5 


121.8 


0.88 


58.9 


<-5 


g--- 




82.6 


0.9 




476 


063200.65+045241.2 


98.002718 


4.878124 


0.1 


3.8 


185.9 


14.2 


1.1 


27.3 


0.89 


12.6 


<-5 




a 


90.6 


1.2 




615 


063209.32+044924.4 


98.038846 


4.823469 


0.1 


7.7 


1685.2 


41.6 


4.8 


17.6 


0.69 


40.0 


<-5 






83.9 


1.0 


CO 


630 


063210.47+045759.6 


98.043639 


4.966572 


0.1 


4.0 


332.6 


18.8 


1.4 


0.3 


0.89 


17.2 


<-5 




a 


88.5 


0.9 


Ol 


743 


063220.77+045303.5 


98.086542 


4.884328 


0.3 


7.0 


25.5 


5.7 


1.5 


2.2 


0.48 


4.1 


<-5 




b 


85.0 


0.9 



Note. — Tabic [2] is published in its entirety in the eleetronic edition of the Astrophysical Journal. The first five sources and the interesting sources shown in Figure 1 are listed here for 
guidance regarding its form and content and for the convenience of the reader. 



Note. — Column 1: X-ray catalog sequence number, sorted by RA. Column 2: lAU designation. Columns 3,4: Right ascension and declination for epoch J2000.0. Column 5: Estimated 



random component of position error, Icr, computed as 



standard deviation of PSF inside extraction region 



Column 6: Off-axis angle. Columns 7,8: Estimated net counts extracted in the 



\/if^ of counts extracted 

total energy band (0.5—8 keV); average of the upper and lower Ict errors on column 7, Column 9: Background counts extracted (total band). Column 10: Estimated net counts extracted 
in the hard energy band (2—8 keV). Column 11: Fraction of the PSF (at 1.497 keV) enclosed within the extraction region. Note that a reduced PSF fraction (significantly below 90%) may 



indicate that the source is in a crowded region. Column 12: Photometric significance computed as 



net counts 



Column 13: Log probability that extracted counts (total 



upper error on net counts ' 

band) are solely from background. Some sources have P^ values above the 1% threshold that defines the catalog because local background estimates can rise during the final extraction 
iteration after sources are removed from the catalog. Column 14: Source anomalies: g — fractional time that the source was on a detector (FRACEXPO from mkarf) is < 0.9 ; e — source 
on field edge; p — source piled up; s ^ source on readout streak. Column 15: Variability characterization based on K-S statistic (total band): a — no evidence for variability (0.05 < Pks)' 
b ^ possibly variable (0.005 < Pk S < 0.05); c ^ definitely variable (P^ S < 0.005). No value is reported for sources with fewer than 4 counts or for sources in chip gaps or on field edges. 
Column 16: Effective exposure time: approximate time the soLircc would have to be observed on axis to obtain the reported number of counts. Column 17: Background-corrected median 
photon energy (total band). 



Table 3 
Chandra Secondary Catalog: Tentative Source Properties 



CO 





Source 




Position 








Ext 


[■acted Counts 








Chai-ac 


tcristics 






Scq 


CXOU J 


Oi ,12000 


*,I2000 


Err 


e 


Net 


ANet 


Bkgd 


Not 


PSF 


Signif 


log Pb 


Anom 


Var 


EffExp 


B„edi„„ 


# 




(dcg) 


(dog) 


(") 


{') 


Full 


Full 


Full 


Hard 


Frac 










(ks) 


(keV) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


(17) 


19 


063124.47+045306.2 


97.851975 


4.885077 


0.8 


8.4 


9.0 


4.9 


10.0 


0.1 


0.90 


1.6 


-2.1 






a 


82.8 


1.5 


28 


063126.33+045728.8 


97.859728 


4.958023 


0.8 


7.4 


8.0 


4.4 


7.0 


3.0 


0.90 


1.6 


-2.2 






a 


84.7 


1.1 


32 


063127.43+045036.4 


97.864312 


4.843457 


0.8 


9.0 


10.3 


5.5 


13.7 


2.8 


0.90 


1.7 


-2.1 






a 


83.0 


1.5 


36 


063128.13+045008.0 


97.867223 


4.835559 


0.8 


9.2 


11.7 


5.8 


15.3 


0.0 


0.90 


1.8 


-2.4 






a 


83.7 


1.3 


52 


063130.85+044847.3 


97.878575 


4.813144 


0.9 


9.7 


13.4 


6.8 


23.6 


0.1 


0.89 


1.8 


-2.2 






a 


83.3 


1.4 


61 


063132.25+050439.0 


97.884410 


5.077512 


0.9 


10.2 


12.3 


6.8 


22.7 


6.9 


0.90 


1.7 


-2.0 






a 


63.1 


2.7 


68 


063133.57+045233.3 


97.889886 


4.875920 


0.6 


6.7 


8.2 


4.1 


4.8 


2.9 


0.89 


1.7 


-2.9 






a 


86.1 


1.8 


71 


063134.14+044958.5 


97.892258 


4.832939 


0.7 


8.3 


10.6 


5.3 


11.4 


0.7 


0.90 


1.8 


-2.5 






a 


84.7 


1.1 


83 


063135.41+045813.0 


97.897576 


4.970293 


0.6 


5.4 


5.8 


3.4 


2.2 


4.5 


0.90 


1.5 


-2.7 






a 


88.9 


4.9 


85 


063135.75+050259.6 


97.898985 


5.049898 


0.7 


8.4 


10.7 


5.3 


11.3 


2.5 


0.90 


1.8 


-2.5 






b 


79.9 


1.6 



Note. — Tablc[3]is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown here for guidance regarding its form and content. 



Estimated random component of position error, la, computed as 



Column 6: Off-axis angle. Columns 7,8: Estimated net counts 



Note. — Column 1: X-ray catalog sequence number, sorted by RA. Column 2: lAU designation. Columns 3,4: Right ascension and declination for epoch J2000.0. Column 5: 

standard deviation of PSF inside extraction region 
\/# of counts extracted 

extracted in the total energy band (0.5—8 keV); average of the upper and lower Icr errors on column 7. Column 9: Background counts extracted (total band). Column 10: Estimated net 
counts extracted in the hard energy band (2—8 kcV). Column 11: Fraction of the PSF (at 1.497 keV) enclosed within the extraction region. Note that a reduced PSF fraction (significantly 



below 90%) may indicate that the source is in a crowded region. Column 12: Photometric significance computed as 



net counts 



Column 13: Log probability that 



upper error on net counts ' 

extracted counts (total band) are solely from background. Some sources have P^ values above the 1% threshold that defines the catalog because local background estimates can rise during 
the final extraction iteration after sources are removed from the catalog. Column 14: Source anomalies: g — fractional time that the source was on a detector (FRACEXPO from mkarf) 
is < 0.9 ; e — source on field edge; p ^ source piled up; s — source on readout streak. Column 15: Variability characterization based on K-S statistic (total band): a — no evidence for 
variability (0.05 < Pks)'i ^ — possibly variable (0.005 < PkS *^ 0.05); c ^ definitely variable {PkS < 0.005). No value is reported for sources with fewer than 4 counts or for sources in 
chip gaps or on field edges. Column 16: Effective exposure time: approximate time the source would have to be observed on axis to obtain the reported number of counts. Column 17: 
Background-corrected median photon energy (total band). 



Table 4 
X-RAY Spectroscopy for Photometrically Selected Sources: Thermal Plasma Fits 



CO 





Source^ 








Spectral Fit'' 






X- 


ray Luminosities'^ 




Notes^ 


Seq 


CXOU J 


Net 


Signif 


log Nh 


kT 


log EM 


log-Ls 


logih 


log-Lh,c 


log Lt 


log it, c 




# 




Counts 




(cm-2) 


(keV) 


(cm-») 






(ergs s-l) 








(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


1 


063114.36+045303.0 


23.3 


3.7 


20.0 


2.0 


53.4 


30.20 


29.96 


29.96 


30.40 


30.40 




2 


063117.04+045228.3 


32.0 


4.5 


21.1 +0.5 


2.0 


-0.2 53.7 +0.2 


30.31 


30.19 


30.20 


30.56 


30.64 




3 


063118.11+045511.8 


12.4 


2.5 


21.6 +1.3 


-2.6 2.8 


53.4 


29.78 


30.03 


30.06 


30.23 


30.39 




4 


063118.29+045223.1 


14.6 


2.7 


21.1 


2.0 


53.3 


29.96 


29.85 


29.86 


30.21 


30.30 




5 


063118.76+045207.6 


9.7 


2.1 


20.8 +1.1 


0.6 


52.9 +0.8 


29.73 


28.33 


28.34 


29.75 


29.85 




6 


063119.27+045110.7 


22.6 


3.6 


20.1 


1.9 


53.4 


30.17 


29.91 


29.92 


30.36 


30.37 




7 


063120.48+045023.4 


98.9 


8.5 


21.5 +0.1 


-0.4 1.7 +0.8 


-0.17 53.8 +0.07 


30.25 


30.18 


30.20 


30.52 


30.73 




8 


063120.89+045003.8 


599.7 


22.8 


21.4 


0.6 


54.4 


31.11 


29.86 


29.89 


31.13 


31.51 


H; HD 46056 


9 


063120.95+045322.7 


14.5 


2.6 


21.0 


3.4 


53.0 


29.70 


29.81 


29.81 


30.06 


30.11 




11 


063121.46+045405.3 


15.6 


2.7 


21.4 


2.0 


53.2 


29.71 


29.70 


29.72 


30.01 


30.16 





Note. — Table |4]is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown here for guidance regarding its form and content. 

''For convenience columns 1—4 reproduce the source identification, net counts, and photometric significance data from Table [2] 

''All fits used the "wabs(apec)" model in XSPEC and assumed 0.3^0 abundances jlmanishi et al.ll200ll : iFeigelson et al.ll2002t ). Columns 5 and 6 present the best-fit values 
for the column density and plasma temperature parameters. Column 7 presents the emission measure for the model spectrum, assuming a distance of 1.4 kpc. Quantities in 
italics were frozen in the fit. Uncertainties represent 90% confidence intervals. More significant digits are used for uncertainties < 0.1 in order to avoid large rounding errors; 
for consistency, the same number of significant digits is used for both lower and upper uncertainties. Uncertainties are missing when XSPEC was unable to compute them or 
when their values were so large that the parameter is effectively unconstrained. Fits lacking uncertainties should be considered to be merely a spline to the data to obtain rough 
estimates of luminosities; actual parameter values are unreliable. 

"^X-ray luminosities are presented in columns 8—12: s = soft band (0.5-2 keV); h = hard band (2-8 keV); t = total band (0.5—8 keV). Absorption-corrected luminosities are 
subscripted with a c; they are omitted when logNjj > 22.5 since the soft band emission is essentially unmeasurable. 

''2T means a two-temperature model was used. H means the fit was performed by hand, usually because the automated fit yielded non-physical results. Names of well-known 
OB counterparts arc listed here for the convenience of the reader. 



Table 5 
X-RAY Spectroscopy for Photometrically Selected Sources: Power Law Fits 







Source^ 








Spectral Fit'" 










X-ray Fluxes'^ 




Notes 




Seq 


CXOUJ 


Net 


Signif 


log Nh 




r 


logJVr 


log is 


logife 


logi/i,c 


logLt 


log -^t,c 






# 




Counts 




(cm-2) 












(photons cm~^ 


s-i) 








(1) 


(2) 


(3) 


(4) 


(5) 




(6) 




(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 




10 


063121.28+045023.8 


50.6 


5.7 


21.5 


-0.9 


1.4 


-0.5 


-5.6 


29.80 


30.45 


30.46 


30.54 


30.62 








15 


063123.59+045141.7 


60.6 


6.3 


21.5 +0.5 


-0.6 


1.1 +0.9 


-0.4 


-5.9 +0.5 


29.63 


30.41 


30.42 


30.47 


30.53 








26 


063125.51+045252.2 


45.5 


5.4 


20.0 




2.0 




-5.8 


29.86 


29.73 


29.73 


30.10 


30.11 








31 


063127.01+050205.0 


14.3 


2.3 


22.3 




1.5 




-6.0 


28.79 


29.92 


29.99 


29.95 


30.17 








84 


063135.68+045322.2 


26.9 


4.1 


-0.7 22.1 +0.4 


-1.2 


2.0 




-5.7 


29.21 


29.98 


30.03 


30.04 


30.33 








89 


063136.33+045251.1 


46.2 


5.7 


-0.6 21.9 +0.3 


-0.9 


1.8 




-5.6 


29.46 


30.15 


30.19 


30.23 


30.43 








91 


063136.42+045602.8 


12.2 


2.5 


22.3 




1.5 




-6.1 


28.70 


29.84 


29.91 


29.87 


30.10 






Co 


94 


063136.49+045959.0 


10.6 


2.2 


21.1 +1.0 




0.8 




-6.8 


28.90 


29.72 


29.72 


29.78 


29.80 






00 


118 


063139.24+050244.7 


19.5 


3.1 


22.1 +0.3 




1.3 




-6.0 


29.02 


30.08 


30.13 


30.12 


30.27 








126 


063139.92+050059.4 


47.1 


5.8 


21.7 +0.5 


-0.8 


1.3 


-0.4 


-5.8 


29.54 


30.31 


30.33 


30.38 


30.47 







Note. — Table[5]is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown here for guidance regarding its form and content. 

''For convenience columns 1—4 reproduce the source identification, net counts, and photometric significance data from Table [2] 

''AH fits used the "wabs(powerlaw)" model in XSPEC. Columns 5 and 6 present the best-fit values for the column density and power law photon index parameters. 
Column 7 presents the power law normalization for the model spectrum. Quantities in italics were frozen in the fit. Uncertainties represent 90% confidence intervals. 
More significant digits are used for uncertainties < 0.1 in order to avoid large rounding errors; for consistency, the same number of significant digits is used for both 
lower and upper uncertainties. Uncertainties are missing when XSPEC was unable to compute them or when their values were so large that the parameter is effectively 
unconstrained. Fits lacking uncertainties should be considered to merely be a spline to the data to obtain rough estimates of luminosities; actual parameter values are 
unreliable. 



■^X-ray luminosities are presented in columns 8—12: s = 
luminosities are subscripted with a c; they are omitted when 



soft band (0.5—2 keV); h = hard band (2—8 keV); t = total band (0.5-8 keV). Absorption-corrected 
og Nh > 22.5 since the soft band emission is essentially unmcasurable. 



Table 6 

Stellar Counterparts 





X-ray Source'^ 
















Opt 


cal/Infrarcc 


Photometry 












Seq 


CxOu J 


USnO Bl.O 


M.ID95 


PS02 BC02 


U 


B 


V 


R 


I 


Hq 


2MASS 


FLAMINGOS 


.1 


H 


K 


PhCcFlg 


# 












(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


ID 


ID 


(mag) 


(mag) 


(mag) 




(1) 


(2) 


(3) 


(4) 


(5) (6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


(17) 


(18) 


1 


063114.36+045303.0 


0948-0096209 










19.33 




17.15 


17.37 




06311429+0453032 






14.49 


13.64 


13.45 


AAAOOO 


2 


063117.04+045228.3 


0948-0096228 










20.75 




18.32 


16.59 




06311696+0452281 






14.94 


14.17 


13.85 


AAAOOO 


3 


063118.11+045511.8 




































4 


063118.29+045223.1 


0948-0096235 










19.47 




16.79 


15.67 




06311825+0452234 






14.53 


13.53 


13.22 


AAAOOO 


5 


063118.76+045207.6 


0948-0096241 


503 




14 


13.36 


13.36 


12.90 


12.38 


12.91 


12.57 


06311881+0452089 






11.98 


11.77 


11.71 


AAAOOO 


6 


063119.27+045110.7 


0948-0096246 










19.39 




16.84 


15.52 




06311925+0451121 






14.35 


13.58 


13.29 


AAAOOO 


7 


063120.48+045023.4 






















06312048+0450239 






13.42 


12.70 


12.46 


AAAcOO 


8 


063120.89+045003.8 


0948-0096261 


454 




15 


7.64 


8.37 


8.22 


8.16 


8.09 


8.13 


06312087+0450038 


063117+050522 


7.84 


7.84 


7.82 


AAAOOO 


9 


063120.95+045322.7 


0948-0096264 










19.50 




17.52 


16.43 




06312100+0453238 


063120+045323 


14.93 


14.14 


13.97 


AAAOOO 


10 


063121.28+045023.8 



































Note. — Tablc[^\vith complete notes is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is slro\vn here for guidance regarding its form and content. 

^Columns 1—2 reproduce the sequence number and source identification from Table [21 and Table ^ Columns 3—6 are the catalogs used for counterparts matching. For convenience, [MJD95]^]VIassey, Johnson, . 
Degioia-Eastwood (1995), [BC02] = Berghofer & Christian (2002), [PS02]=Park & Sung (2002). Columns 7-12 give available optical photometry. Columns 13-17 provide NIR identifications and JHK photometry froi 
FLAMINGOS (§3.1). Column 18 lists the 2MASS photometric quality flags (Cutri et al. 2003). In the note on individual sources, [OlSlJ^Ogura & Ishida (1981), [LR04]=Li & Rector (2004), ProbMem=Probability of bein 
a cluster member from proper motion data: [M82]— Marschall et al. (1982), [D06]— Dias et al. (2006). 

^2"from [OI81] 76=[BC021 14 

^=HD 46056=BD+04 1291=MWC 808=[OI811 84 (08V); FUSE spectrum available 



Table 7 
x-ray properties of cataloged ob stars in ngc 2244 





Optica 


/IR Properties 












X-ray Properties 






Name 


SpTy 


2MASS 


K 


logito! 


Seq 


Ac/) 


NetCts 


logiVff 


kT 


logih 


log it, c 








(mag) 


(L©) 


# 


{") 




(em-2) 


(keV) 


(ergs~l) 


(ergs~l) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


HD 46223 


04V((f)) 


06320931+0449246 


6.68 


5.7 


615 


0.2 


1685 


21.6 


0.3 


29.71 


32.38 


HD 46150 


05V((f)) 


06315551+0456343 


6.44 


5.5 


373 


0.1 


3589 


21.4 


0.6+0.2 


30.48 


32.34 


HD 46485 


07V 


06335094+0431316 


7.45 


5.2 


RMC 164 


0.3 


310 


21.6 


0.3+0.9 


30.46 


32.05 


HD 46056 


08V((f)) 


06312087+0450038 


7.82 


5.0 


8 


0.3 


600 


21.4 


0.6 


29.86 


31.51 


HD 46149 


08.5V((f)) 


06315253+0501591 


7.25 


4.9 


311 


0.1 


299 


20.0 


0.7 


29.62 


31.04 


HD 258691 


09V((f)) 


06303331+0441276 


7.93 


4.8 


NFOV 














HD 46202 


09V((f)) 


06321047+0457597 


7.72 


4.8 


630 


0.1 


333 


21.2 


0.3 


28.75 


31.11 


HD 259238 


BOV 


06321821+0503216 


10.28 


4.5 


727 


1.0 


36 


21.8 


0.5 


28.74 


30.65 


HD 46106 


B0.2V 


06313839+0501363 


7.62 


4.4 


107 


0.2 


186 


21.4 


0.4+1.4 


29.81 


30.94 


MJD95 


B0.5V 


06313708+0445537 


12.20 


4.3 


NFOV 














HD 259135 


B0.5V 


06320061+0452410 


8.12 


4.3 


476 


0.7 


186 


20.6 


1.6 


30.04 


30.57 


IRAS 06309+0450 


B0.5V 


06333749+0448470 


8.64 


4.3 


NFOV 














HD 259012 


BIV 


06313346+0450396 


8.79 


4.0 


66 


0.6 


275 


21.1 


2.5 


30.75 


31.12 


HD 259105 


BIV 


06315200+0455573 


8.95 


4.0 












<28.3 


<28.7 


BD+04° 1299s 


BlIII 


06320613+0452153 


9.38 


4.0 


... a 










<28.3 


<28.7 


HD 46484 


BIV 


06335441+0439446 


6.86 


4.0 


NFOV 














BD+05°1281B 


B1.5V 


06315893+0455398 


9.74 


3.7 


448 


0.1 


16 


21.1 


2.2 


29.58 


30.01 


HD 259172 


B2V 


06320259+0505086 


10.08 


3.5 


NFOV 














OI81 345 


B2 


06330656+0506034 


11.20 


3.5 


NFOV 














BD+04°1295p 


B2.5V 


06313146+0450596 


10.24 


3.4 


53 


0.5 


28 


21.5 


2.0 


29.57 


30.03 


OI81 130 


B2.5V 


06314789+0454181 


10.87 


3.4 












<28.3 


<28.7 


OI81 190 


B2.5VI1 


06315891+0456162 


10.54 


3.4 












<28.3 


<28.7 


OI81 172 


B2.5V 


06320984+0502134 


10.43 


3.4 












<28.3 


<28.7 


OI81 274 


B2.5V 


06322424+0447037 


10.58 


3.4 












<28.3 


<28.7 


OI81 392 


B2.5V 


06335056+0501376 


9.99 


3.4 


NFOV 














OI81 194 


B3 


06321548+0455203 


10.96 


3.2 


697 


0.1 


62 


21.2 


3.6 


30.16 


30.45 


MJD95 


B3V 


06322249+0455342 


13.48 


3.2 












<28.3 


<28.7 


HD 259268 


B3 


06322304+0502457 


10.33 


3.2 












<28.4 


<28.8 


HD 259300 


B3Vp 


06322939+0456560 


9.34 


3.2 


815 


0.3 


60 


21.5 


1.4 


29.88 


30.54 


OI81 334 


B3 


06325179+0447161 


11.53 
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Note. — Column 1: This list is obtained from Appendix A of TFM03 which gives optical cross-identifications, positions, and spectral types. The stars are listed 
first in order of decreasing mass, and then by right ascension. OI— [Qgura &: Ishida i 19811 ). Column 2: Spectral types are from Ogura &; Ishida (1981) and Massey 
et al.(1995). 

Colunnns 3—4: Source numbers and i^'-band magnitudes are from the 2MASS All-sky Point Sources Catalog . 

Column 5: Bolometric luminosities are estimated from calibrations of Lj^^i with spectral type. lMartins et al.| J 2 0051) for 03— 09.5 stars and lde Jager Sz Nieuwenhuijzenl 

il987|) for B stars. No use is made of available photometry. 

Column 6: Chandra source number, from Table [^ NFOV— object is not covered in the FOV. ...^non detection. HD 46485 is observed in Rosette Field 4 (RMC 

source #164, Paper II). 

Column 7: Offset (in arcscconds) between the Chandra and 2MASS sources. 

Columns 8—12: X-ray properties from Tablc[^ extracted counts after background subtraction; column density and plasma energy from fits to the ACIS spectra (: 

in log A'"/f values arc approximated from median energy ( Fcigelson et al. 2005^; : in kT are assumed); observed hard band luminosity (2 — 8 koV); inferred total band 

luminosity corrected for absorption (0.5 — 8 keV). Upper limits to luminosities are estimated using the faintest sources in Tablc[^and scaled with the corresponding 

exposure time. 

^BD-(-04° 1299s was reported as a detection in TFM03 (20ks observation). However the deep observation resolved this source into two X-ray sources. Both are 
separated by --^2 from the optical position. Therefore we do not report BD + 04° 1299s as a detection here. 
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Fig. 1. — (a): A large scale ( ^ 2° x 1.5°) view of the Rosette star forming complex in CO J =1 — 



emission from iHeyer et al.l (120061 ). The multiple ACTS FOVs (polygons) and the extent of the 
NGC 2244 clus t er (c ircle) are shown. Squares mark the embedded clusters in the RMC with 
Phelps &; Ladal (|l997l ) sequence numbers, (b): A 45' x 45' DSS2 i?-band image of the Rosette 
Nebula. Ah known O stars in the FOV that belong to the NGC 2244 cluster (TFM03 Table 6) 
are labeled by circles. Four stars are marked as squares: V578 Mon is an eclipsing binary; HH 1 
is a stellar microjet; ACIS #919 is a candidate massive star; the visually brightest star HD 46241 
(KOV) is foreground. These objects are described further in the text. 
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Fig. 2.— (a): A merged 94 ks ACIS-I image of the NGC 2244 cluster from ObsIDs 1874 and 3750 
(outlined by two 17' x 17' boxes) with reduced resolution (binned by 2 pixels). The two ObsIDs 
have slightly different roll angles, (b): X-ray composite image created from csmooth for the merged 
fields. Blue intensity is scaled to the soft (0.5-2 keV) X-ray emission, green intensity is scaled to the 
hard (2-7 keV) X-ray emission, (c): Same as (b) but the scaling emphasizes soft diffuse emission 
and red intensity is scaled to the DSS i?-band optical emission. 
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Fig. 2. — Continued. 



44 



^Fi 







\ 

J- 






.■■ ■. ■■■■.®. 



^ 



,. tj^ W »; o 



V ^> 




^© 



^ ^^^ 




Fig. 3. — 30" X 30" ACIS unbinned image (top left), reconstructed image (top right), the Ha image 
(bottom left), and the 2MASS-i^s image (bottom right) of the central region around the 05V 
star HD 46150. Several new sources are resolved within 5" of the dominant star by the Chandra 
observation. 
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Fig. 4. — Lightcurves of sources with more than 500 counts that are significantly variable (Pks ^ 
0.005). The ACIS sequence numbers and binsizes are marked. 
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Fig. 5. — NIR J — H vs. H — K color-color diagram for 617 Chandra stars with high-quality 
photometry, from combined FLAMINGOS and 2MASS data (error in both J — H and H — K colors 
< 0.1 mag). The (light and dark) green circles and red triangles represent sources with significant 
X-band excess and sources without excess, respectively. The dark green circles represent three 
Class I objects and are labeled with their sequence numbers from Table 2. The five blue triangles 
are foreground stars. Stars using 2MASS photometry are indicated with black circles. The black 
solid and long-dash lines denote the loci of MS stars and giants, respectively, from Bessell & Brett 
(1988). The purple dash dotted line is the locus for classical T Tauri stars from Meyer et al. (1997), 
and the cyan solid line is the locus for HAeBe stars from Lada & Adams (1992). The blue dashed 
lines represent the standard reddening vector with crosses marking every Ay = 5 mag. Most 
Chandra sources are located in the reddening band defined by the left two dashed lines associated 
with Class III objects (triangles). To the right of this reddened band are 38 IR-excess sources. 
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Fig. 6. — NIR J vs. J—H color-magnitude diagram using the same sample and symbols as FigureO 
except that known O stars are denoted as the red filled triangles. ACIS source n umbers are r narke d 
for some stars. The purple solid line is the 2 Myr isochrone for PMS stars from lSiess et al.l ( 200CI ). 
The dash dotted line marks the location of Zero Age Main Sequence (ZAMS) stars. The blue 
dashed lines represent the standard reddening vector with asterisks marking every Ay = 5 mag 
and the corresponding stellar masses are marked. 
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Fig. 7. — The cumulative distribution of the source hardness indicator, median photon energy, 
for Chandra sources with identified ONIR counterparts (solid line) and those without counterparts 
(dashed line). The sources with identified ONIR counterparts are considerably softer than the 
latter group. 
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Fig. 8. — The spatial distribution of Chandra sources without identified ONIR counterparts. The 
background is the DSS i?-band image. ACIS-I FOVs are shown. 
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Fig. 9. — (a): X-ray luminosity function (XLF) constructed from the absorption corrected hard 
band (2.0-8.0 keV) X-ray luminosity L^^c for the unobscured NGC 2 244 population (solid line), 
and the COUP ONC unobscured cool stars population (dashed line, iFeigelson et al.ll2005l ). The 
vertical line denotes the estimated completeness limit for the NGC 2244 population, (b): The 
cumulative distribution of X-ray luminosities for the unobscured N GC 2244 X-ray sourc es as well as 
the unobscur ed populations in C OUP. Cep B/0B3b (age ^l-3 Mvr.lCetman et ahl boOfil'l . NGC 6357 
(age~l Mvr. IWang et al.ll2007l ). and M17 (age~l Mvr. iBroos et al.ll2007l ). The distributions for 
NGC 2244, NGC 6357, Cep B, and M17 are truncated at their completeness limits, log Lh,c ~ 29.2, 
30.4, 29.3, and 30.4 ergs s~^, respectively. 
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Fig. 10.— Comparison between IMFs of the NGC 2244 X-ray stars (solid line) and COUP ONC 
stars (dashed line) using NIR photometry derived stellar masses. The dotted lines show the ONC 
IMF (JMuench et al.ll2002l ) and its scaled version to match the NGC 2244 IMF. The dot-dashed hne 
is the NGC 2244 IMF estimated using the KLF derived from 2MASS data. The arrows indicate 
the approximate mass completeness limits for the NGC 2244 and ONC X-ray stars, and the mass 
completeness limit for the KLF derived from 2MASS data. Note that the 2MASS completeness 
limit is the same as our Chandra completeness limit in stellar mass. The incomplete bins in the 
NGC 2244 KLF/IMF are omitted. 
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Fig. 11. — The stellar surface density map for the unobscured stellar population in NGC 2244. The 
cluster shows a spherical structure that extends 8 arcmin in diameter. The highest concentration 
of stars is around R.A.=06^31™55^ Dec.=+04°56'34", the location of HD 46150. A secondary 
density enhancement is seen centered at R.A.=06^31™56^ Dec.=04°54'10". The 04 star HD 46223 
is mostly isolated. 
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Fig. 12.— The KPNO Ha image (JLi fc Rectoil l2004l ) and 2MASS Ks image of the secondary 
overdensity of stars seen in Figure [TTJ The overlaid circles mark the ACIS sources. The 3.5' x 3.5' 
images are both centered at RA=06^31'^56',Dec.=04°54'l6". 
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Fig. 13. — The observed radi al density profiles of the NGC 2244 cluster, the ONC from COUP 
studies (IFeigelsoii et al.ll2005l ). and the NGC 6357 region from our Chandra/ ACIS observation 
( Wang et al.1 20071). The histogr am shows the radial density profile of the ONC from ONIR studies 
Hillenbrand Sz Hartmannlll998l ). la Poisson error is shown for the COUP ONC. 



55 







1 








' 


1 


1 1 1 1 1 1 


1 






- 






'■, 






Ntc 6:is-? '.'.' 


- 






_ 






' \ 


■ 




t-\f-C '■■'■•A *■ 


_ 






---. \ 


t 






"k- 




COUP ■ - 








■\ 


\ 






H 












^ 




\ 




-M.-: 








y 




\ 

^ 




\ 






r 








tOQO 




'•. 




\ 




»■ -.■ 




1 


:1 




- 








\ 


\ 




- 


f 




_ 








\ 


'^- 




_ 


M 




. 














^ 


a 




■ — 










-^SSt^\ 






•^ 










m 














^^^"-O^ 


* 




,i_ 














■■^ ^ 


lV 




■-«; 














\ 


^l^v.^ 




a 


too 


— 










V 


^^^\rf- 


— 


■si; 




_ 










"N- 


>vV 


^ 


W: 


10 


- 1 













1 1 ll '^l'-- 1 1 


r 



U1 



m.o. 



Rfidiai 'dsstonte (pc) 



Fig. 14. — The radial density profiles (thick lines) in physical scales for the same three clusters, 
where the stellar densities of NGC 2244 and NGC 6357 have been scaled to 1.2 times and 5 times 
the ONC population, respectively, based on the XLF analysis. The thin line represents a King 
model profile for the outer portion of NGC 2244. 
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Fig. 15. — (a): Spatial distribution of the massive stars (NIR estimated mass M > SMq, filled stars) 
and the low mass stars (M < SMq, open diamonds) in NGC 2244, using our X-ray-selected sample, 
(b): The cumulative radial distributions for the massive stars (NIR estimated mass M > 8Mq, 
dashed line) and the low mass stars (M < 2Mq, solid line). 
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Fig. 16. — The spatial distribution of 38 sources with significant NIR color-excess (circles). The 
northern part of the nebula seems deficient in NIR excess sources. The boxes outline the FOVs of 
the multiple ObsIDs. 
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Fig. 17. — Spectral fits to X-ray spectra of six O and early B stars. Source name, source counts, 
and fit parameters are marked in each panel. The two model components are shown as the dotted 
lines for spectra that are best fit with two-temperature thermal plasma models. 
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Fig. 18. — The Lx vs. Lboi relation for X-ray detected O and early B stars. T he samples are from 
the NGC 2244 cluster (filled triangles; this work, Tabl e [7D, the ONC (crosses, [Stelzer_el_alJ[2005|), 
and the ma,ssive star forming regions M17 (diamonds. iBroos et al.l 120071 ) and NGC 6357 (pluses, 



fromlBroos et al. 



Wang et al.l 12007 ). Upp er limits are marked as arrows. The bolometric luminosities are adopted 



([20071) to be consistent. 
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Fig. 19. — (a): H^ image of the neighborhood of ACIS source #919, near a molecular pillar, (b): 
The X-ray image of the same region. The 2.5' x 2.5' images are both centered at the X-ray bright 
star #919. (c): The spectral fit to the X-ray spectrum of source #919 with logNn = 21.2 and a 
hard {kT = 3.3 keV) plasma. 
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Fig. 20. — X-ray light curve of the echpsing binary V578 Mon (lower panel), together with the 
optical light curve (6-band photometry with Icr error-bars taken from Hensberge et al. 2000) in 
or bital phases (upper pan el). The orbital phase is calculated using the eclipse ephemeris reported 
in iHensberge et al.l (J200d ). 



